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Microwave Radiometer Studies of 
Atmospheric Water Over the Oceans 

A Technical Report in Two Volumes 


P.I.: Kristina B. Katsaros 
Department of Atmospheric Sciences 
UNIVERSITY OF WASHINGTON 
Seattle, WA 98195 
(206) 543-1203 
(206) 543-0308 


Introduction 

Since the Seasat satellite carried the Scanning Multichannel Microwave Radiometer 
(SMMR) into space in July of 1978, shortly followed by the SMMR on Nimbus 7, which 
operated for almost a decade, a new type of data source on atmospheric water vapor and 
other meteorological parameters has been available for analysis of weather systems over the 
ocean. Since 1987 we have had the Scanning Multichannel Microwave/Imager (SSM/I) 
instrument on Defense Meteorological Satellites providing similar data. In this two volume 
technical report we present a collection of our work performed over the last three years 
under NASA grant NAGW-1688 administered by Dr. James Dodge of the Mesoscale 
Program at NASA Headquarters, and NASA grant NAG5-943 administered by Dr. 
Ramesh Kakar of the Global Processes Program. Some of the chapters are journal articles, 
others are conference proceedings and abstracts. A major contribution from our group is 
the Ph.D. thesis of Grant W. Petty, 1990, produced separately as Final Technical Report 
on NASA grant NAG5-943. It is available upon request. 

Volume I of this report contains work mainly supported by NASA grant NAGW-1688. 
The first is an article for publication in Monthly Weather Review: 

Miller, D.K. and K.B. Katsaros, 1992: Satellite derived surface latent heat 
fluxes in rapidly intensifying marine cyclones. Mon. Wea. Rev., 
(scheduled publication: June). 
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It is derived from the following Masters of Science thesis: 

Miller, Douglas K., 1990: Estimation of Surface Latent Heat Flux Patterns 
Preceding a Rapidly Intensifying Cyclone Derived from the Special Sensor 
Microwave! Imager. M. Sci. Thesis, Department of Atmospheric Sciences, 
University of Washington, Seattle, WA 98195. [Not included in this 
report.] 

In this thesis a method for evaluating the preconditioning of the marine environment via 
surface moisture flux 24 hours in advance of rapid deepening of a cyclonic storm is 
explored with data from a 1988 pre-ERICA (Extratropical Rapidly Intensifying Cyclones 
over the Atlantic) storm. We plan to carry on this project with studies of the 1989-90 
ERICA storms (Hadlock and Kreitzberg, 1988). 


One of our goals was to make full use of the new tool for studying marine weather systems 
that the SSM/I represents by applying our algorithms and versatile graphics presentation to 
the task at hand. The hardware and software system were developed by G.W. Petty while 
a graduate student with funds supplied by a special NASA grant (NAS8-36473, Cyclonic 
Storms). 

Collaboration with the following colleagues and visiting student researchers have enabled 
us to advance along several avenues at once: Dr. Nelly Mognard at the University of 

Puget Sound; Dr. Chantal Claud of Ecole Polytectinique Paliseau, student of Professor 
Alain Chedin; and Ms. Aurelie Sand of the University of Paris VII, student of Professor 
Claude Klapisz. Results of the work presented at the European Geophysical Society 
meeting in Copenhagen, April 23-27, 1990, were reported on in our Third Annual Report, 
1990. The following new work was presented at the Royal Meteorological 
Society/American Meteorological Society joint meeting in London, September 3-7, 1990: 


Claud, C., N. Scott, A. Chedin, K. Katsaros and G. Petty, 1990: Mesoscale 
Meteorology in Polar Regions from a Coupling of NOAA (TOVS) and 
DMSP (SSM/I) Data. In Preprint Volume, Fifth Conference on Satellite 
Meteorology and Oceanography, September 3-7, 1990, London, U.K. 
Published by the Am. Meteorol. Soc., Boston, MA. 

Mognard, N.M. and K.B. Katsaros, 1990: Evolution of Atmospheric Fronts 
over the Ocean as Observed with the Special Sensor Microwave/Imager and 
the Geosat Altimeter. In Preprint Volume, Fifth Conference of the Royal 
Meteorological Society and the American Meteorological Society, London, 
U.K., September 3-7, 1990. 
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The following article which will appear in Tellus is an extension of the London presentation 
by Claud et al. (1990): 


Claud, C., K.B. Katsaros, G.W. Petty, A. Chedin and N.A. Scott, 1992: A 
cold air outbreak over the Norwegian Sea observed with the Tiros-N 
Operational Vertical Sounder (TOVS) and the Special Sensor 
Microwave/Imager (SSM/I). Tellus, (in press). 

Follow-up papers to the conference presentations by Mognard, and by Sand, are also being 
prepared. 


Earlier work supported by the Mesoscale Program of NASA appeared this year as a journal 
article: 


McMurdie, L.A. and K.B. Katsaros, 1991: Satellite-derived precipitable water 
distribution in oceanic midlatitude storms: variation with region and season. 
Mon. Wea. Rev., 119, 589-605. 


It is also included in this report, Vol. I. Collaborative work with Scandinavian colleagues 
in relating the SMMR or SSM/I water parameters to similar quantities extracted from 
numerical models of the atmosphere have resulted in a University of Stockholm Technical 
Report: 


Katsaros, K.B., G.W. Petty and U. Hammarstrand, 1990: Atmospheric Water 
Parameters in Mid-Latitude Cyclones Observed by Microwave Radiometry 
and Compared to Model Calculations. Technical Report, Meteorological 
Institute, University of Stockholm (in press). [Not included in this 
collection but available on request] 


An extension of ideas for this report were carried out with a later version of the Stockholm- 
Bergen mesoscale atmospheric model. A short article for Tellus is included in this report, 
Vol. I: 


Raustein, E„ H. Sundqvist and K.B. Katsaros, 1991: Quantitative comparison 
between simulated cloudiness and clouds objectively derived from satellite 
data. Tellus, 43A, 306-320. 

Volume II of this technical report contains work mainly supported by NASA grant NAG5- 
943 and its predecessor NASA grant NAG5-354. We recognized early on the importance 
of microwave radiometer estimates of precipitation over the oceans, but also the difficulty 
in evaluating the accuracy of such estimates since few comparison data exist. Radar data 
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provide a source for comparison, but the difference in resolution between the satellite 
estimates of the order of 35 km and the radar estimates of the order of 0.5 km are a 
problem. Two articles concerned with these problems have been produced and are 
reproduced here: 

Petty, G.W. and K.B. Katsaros, 1990: Precipitation observed over the South 
China Sea by the Nimbus 7 Scanning Multichannel Microwave Radiometer 
during Winter MONEX. J. Appl. Meteor., 29, 273-287. 

Petty, G.W. and K.B. Katsaros, 1992: Nimbus 7 SMMR precipitation 
observations calibrated against surface radar during TAMEX. J. Appl. 

Met., (scheduled publication: May). 

The research performed by Petty was also presented at the American Meteorological 
Society Satellite Conferences in San Diego, 1989, and London, 1990. The London 
conference proceedings are: 


Petty, G.W. and K.B. Katsaros, 1990: New geophysical algorithms for the 
Special Sensor Microwave/Imager. Preprint Volume, Fifth International 
Conference on Satellite Meteorology and Oceanography, September 3-7, 
1990, London, England. [Included in this report.] 

Petty, G.W. and K.B. Katsaros, 1990: Nimbus 7 SMMR precipitation 
observations calibrated against surface radar during TAMEX. J. Applied 
Met., (submitted). Also appearing in abbreviated form in Preprint Volume, 
Fifth International Conference on Satellite Meteorology and Oceanography, 
September 3-7, 1990, London, England. [Not included here; overlaps with 
the article above with the same title.]" - 


The thesis work of Petty will appear as a series of articles in the near future. One has been 
submitted. 

A review of the legacy that the Seasat Satellite provided for all of the work in our group on 
microwave radiometry and other work related to the scatterometer and altimeter instruments 
has appeared in the Bulletin of the American Meteorological Society: 

Katsaros, K.B. and R.A. Brown, 1991: Legacy of the Seasat Mission for 
studies of the atmosphere and air-sea-ice interactions. Bull. Am. Met. Soc., 

72, 967-981. 
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Three papers are accepted for presentation at the Sixth International Conference on Satellite 
Meteorology and Oceanography, in early 1992. 


McMurdie, L. and K.Katsaros, 1992: Satellite Derived Integrated Water Vapor 
and Rain Intensity Patterns: Indicators of Rapid Cyclogenesis. In Preprint 
Volume, 6th Conference on Satellite Meteorology and Oceanography, Jan. 
5-10, 1992, Adanta, GA. 

Mognard, N.M. and K.B. Katsaros, 1992: Comparison of Wind Speed 
Measurements over the Oceans with the Special Sensor Microwave/Imager 
and the Geosat Altimeter. In Preprint Volume, 6th Conference on Satellite 
Meteorology and Oceanography, Jan. 5-10, 1992, Atlanta, GA 

Petty, G.W. and K.B. Katsaros, 1992: Morning-Evening Differences in Global 
and Regional Oceanic Precipitation as Observed by the SSM/I. In Preprint 
Volume, 6th Conference on Satellite Meteorology and Oceanography, Jan. 
5-10, 1992, Atlanta, GA. 


The preprints of these conference papers are the final contribution in Volume n. 


References: 

Abbey, R.F., C.W. Kreitzberg and R. Hadlock, 1987: Experiment on Rapidly 

Intensifying Cyclones over the Atlantic, ERICA, Overview Document, Field 
Research Project, Version 1. March 1987, available from the ERICA Project Office, 
Battelle Memorial Institute, Richland, WA. 

Hadlock, R. and C.W. Kreitzberg, 1988: The Experiment on Rapidly Intensifying 
Cyclones over the Atlantic (ERICA) field study-objectives and plans. Bull. Am. Met. 
Soc., 69, 1309-1320. 

Petty, G.W., 1990: On the Response of the Special Sensor Microwave/Imager to the 
Marine Environment - Implications for Atmospheric Parameter Retrievals. Ph.D. 
Thesis, October 1990, Dept, of Atmospheric Sciences, University of Washington, 
Seattle, WA 98195. NASA Tech Report. 


vi 





Satellite-derived surface latent heat fluxes 
IN A RAPIDLY INTENSIFYING MARINE CYCLONE 


Douglas K. Miller* 

Kristina B, Katsaros 

Department of Atmospheric Sciences AK-40 
UNIVERSITY OF WASHINGTON 
Seattle, WA 98195 
Tel: (206)543-1203 
Fax: (206)543-0308 


Monthly Weather Review 
In press (publication date June 1992) 


* Present affiliation: 

Department of Earth and Atmospheric Sciences 
PURDUE UNIVERSITY 
West Lafayette, IN 47907 






Abstract 


The aim of this article is to estimate surface latent heat fluxes in the vicinity of a 
rapidly deepening cyclone before and during its period of most rapid intensification. 
We do this with a bulk parameterization scheme (Liu, Katsaros, and Businger, 
1979) and remotely sensed input data. 

We investigate a method for estimating the difference in specific humid- 
ity between the surface and a ten meter height using the Special Sensor Mi- 
crowave/Imager (SSM/I)-observed integrated water vapor field and a sea surface 
temperature analysis. 

The surface latent heat flux fields generally have estimated errors below 40% 
south of 40° N and outside the region of high integrated water vapor values as- 
sociated with frontal bands. The method of estimating surface latent heat fluxes 
for the case study was found to be usable in most regions of the Northwest At- 
lantic Ocean except for those locations directly adjacent to coastlines in instances 
of offshore flow and in the vicinity of surface fronts. 
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1. Introduction 

Cyclonic storms that form and move along the east coast of North America 
during the winter months often bring significant amounts of precipitation and 
strong winds to New England and the Maritime Provinces. The need to better 
understand and predict these storms is made clear by the direct effect that they 
have on mankind. In some cases they have inflicted a high cost in life and property. 
Ships have experienced extensive deck damage or been sunk, fishing has been 
disrupted, and an oil rig was once destroyed with eighty- four lives lost. 

Those cyclones which fulfill a specific deepening criterion (i.e. 1 mb hr -1 for 
24 hours, geostrophically adjusted to 60° latitude) are classified as ‘rapid deep- 
oners’ (Sanders and Gyakum, 1980). These storms, also referred to as ‘explosive 
deepeners’, have been the focus of much research in the past decade. A statistical 
analysis of cyclone deepening rates by Roebber (1989) indicates that the Atlantic 
and Pacific cyclones in his sample have a significantly greater variance in cyclone 
deepening rates when compared to the rates for storms located over the continental 
U.S and Canada. Frequency distributions of the 24 hour adjusted maximum deep- 
ening rates produced by Roebber demonstrate that the marine cyclones show a 
greater skewness towards strong deepening than does the continental sample, indi- 
cating that an oceanic cyclone is statistically more likely to have a large deepening 
rate and be classified as a rapid deepener. 

The importance of latent heat release in cyclogenesis is well documented. Reed 
and Albright (1986) examined a case of explosive cyclogenesis in the Eastern Pacific 
and concluded that condensational heating was the only diabatic heating process 
significantly contributing to the cyclonic development. In a numerical simulation 
of the same storm, Kuo and Reed (1988) confirmed that latent heat release and 
its interaction with baroclinity accounted for half of the deepening in the control 
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run. Reed et al. (1988) studied three storms in the North Atlantic and concluded 
that latent heat release was responsible for between forty and fifty percent of 
the deepening in two storms which were well predicted by the European Centre 
for Medium Range Weather Forecasts (ECMWF) T106 operational model. A 
conclusion of the relative importance of latent heat release for the third storm was 
not possible since it was not well predicted. Chen et al. (1983) and Gyakum (1983) 
also demonstrate the important effect of latent heat release as a contribution to 
cyclonic development. Kuo et al. (1991a) simulated the now famous QE-II storm 
and found that latent heat release significantly modified the warm frontal structure 
of the storm to enhance the adiabatic secondary circulation. 

The vertical transport of water vapor at the surface (moisture flux or latent 
heat flux) can affect cyclogenesis indirectly. The upward transport of water vapor 
decreases low and mid-level static stability through the release of latent heat as 
noted by Reed and Albright (1986). Recent work by Kuo and Low-Nam (1990) 
and Kuo et al. (1991b) suggests that surface fluxes (both sensible and latent heat 
fluxes) can actually have a damping effect during the time of rapid intensification 
of explosively deepening storms by decreasing the low level baroclinity since cold 
air tends to be warmed by the surface fluxes and vice versa. Kuo et al. (1991b) 
tested the influence of the surface fluxes by modeling the same storm under varying 
surface flux conditions. Their results indicate that surface energy fluxes during 
the 24 hours that preceded the period of most rapid intensification had a much 
stronger effect on subsequent storm development than did the fluxes that existed 
as the storm rapidly intensified. The fluxes act to precondition the environment 
through moistening and destabilization of air at the surface thereby making water 
vapor available to the developing storm. As the storm spins up, the water vapor 
is converted to cloud water and/or rain and latent heat is released in the process. 
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Latent heat release, particularly in the warm frontal clouds, can enhance the storm 
development by inducing strong low-level convergence and vorticity spin-up as 
observed in the case study presented by Kuo and Reed (1988). 

A preliminary field test for the Experiment on Rapidly Intensifying Cyclones 
over the Atlantic (ERICA, Hadlock and Kreitzberg, 1988) was conducted in late 
January 1988 in order to evaluate the quality of participating observing systems, 
data networks, field operations plans, and data management. Most of the mea- 
surement activity was concentrated in the period 25-26 January and the cyclone 
that developed during this period over the Northwest Atlantic is the focus of this 
investigation. 

Explosively deepening storms are predominantly maritime events as suggested 
by Roebber (1989). The study of cyclones over oceans is made difficult by the 
poor surface observing network currently existing in oceanic regions. Useful data 
is sparse. Instruments mounted on satellites have helped to provide the research 
community with additional data over oceans. For a long time, most satellite in- 
formation has come from images obtained in the visible and infrared wavelengths. 
However, these instruments are limited in the amount of below- cloud- top informa- 
tion they can yield by the thickness of clouds in the field of view. This certainly 
presents serious limitations for studies involving thick cloud formations such as 
those commonly observed in cyclones. Sensors receiving radiation in the microwave 
wavelengths, on the other hand, are capable of giving columnar atmospheric in- 
formation in most cases when significant cloud amounts are present in the field of 
view. The Special Sensor Microwave/ Imager (SSM/I) is one microwave instrument 
currently in operation which can observe oceanic regions. It will be introduced in 
greater detail in section 2. 
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The emphasis. of this study will be on evaluating a method of estimating surface 
latent heat fluxes using data from the SSM/I, sea surface temperature (SST) and 
surface pressure data and to apply the method in order to examine surface latent 
heat flux (LHF) patterns observed in the Western Atlantic both before and during 
the January 1988 case. Since this storm occurred during the pre-ERICA field test, 
a larger set of observations exist than would be expected for most winter storms 
over the ocean. A secondary objective of this study is to develop methodology 
which can be applied to analysis of the ERICA cyclones. SSM/I data for the 
ERICA cyclones was not yet available when this study commenced. The 24-26 
January 1988 case will be referred to hereafter as the ‘pre-ERICA storm’ (Hadlock 
and Kreitzberg, 1988 and Neiman et al., 1990). 

Since the SSM/I instrument doesn’t directly observe moisture flux fields, it 
becomes necessary to use a method by which the moisture flux field is computed 
given a set of mean quantities, the so called bulk parameterization method. Bulk 
parameterization requires that the specific humidity at a ten meter height and at 
the surface (Qa and Qs, respectively) be known. In this study we will develop a 
statistical method for inferring (Qs-Qa) given the SSM/I derived integrated water 
vapor (IWV) fields along with independent sea surface temperature and surface 
pressure fields. This method utilizes the high correlation between the IWV amount 
and the specific humidity at ten meters under certain conditions. The high correla- 
tion is observed because boundary layer air contains 50% or more of the integrated 
water vapor (Liu, 1984). 
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2. Data Sources 

a. The SSM/I Instrument and Its Applications 

In this study, we rely on satellite data derived from the SSM/I. The SSM/I 
is a passive microwave imager which observes electromagnetic radiation at fre- 
quencies that fall between strong absorption bands within the microwave portion 
of the spectrum. As a result, the radiances observed by each SSM/I channel is 
a combination of surface emission or reflection and the emission and attenuation 
occurring through the atmosphere. The SSM/I has seven channels which con- 
sist of two orthogonal polarizations at three frequencies and one polarization at 
a fourth frequency. A geophysical parameter can be theoretically estimated given 
its unique contribution to the total radiant energy received at each SSM/I chan- 
nel. The atmospheric parameters that can be observed by the SSM/I consist of 
near surface wind speed, IWV, cloud liquid water, and an index for rain/no rain 
regions. However, retrieval of all of the listed parameters is not possible at all 
times. Difficulty occurs when the SSM/I is viewing land regions or areas where 
there exist significant amounts of liquid water in the field of view (Petty, 1990). 

The SSM/I IWV and near surface wind speed algorithms used in this study 
were developed and modified by Petty and Katsaros (1990). The near surface wind 
speed algorithm is a modified version of the Goodberlet et al. (1989) algorithm. 
These algorithms were derived using statistical regression of SoM/I brightness 
temperatures (Tg’s) against in situ observations of the respective parameter. The 
algorithms are sensitive to the presence of precipitation and result in less accurate 
parameter estimations as the amount of precipitation in the SSM/I field-of-view 
increases. They become unusable when certain Tb threshold values are exceeded. 
Pixels where the near surface wind speed algorithm Tb threshold values have been 
exceeded are referred to as ‘rain flagged’ pixels. 
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b. Other Data Sources 

By mere coincidence, the SSM/I orbit results in overpasses occurring over the 
North and West Atlantic close to the synoptic reporting times of 0000 and 1200 
UTC. The SSM/I overpasses that occurred in January 1988 to within one hour 
of 0000 or 1200 UTC were chosen for fitting (Qs-Qa) computed from 124 ship of 
opportunity and island station (S/I) surface observations to SSM/I-observed IWV 
and independent SST analyses. All observations that were used to compute the 
statistical fit are located in the North Atlantic from 20 to 50° N and from 55 to 
85° W. The SSM/I overpasses used to compute surface latent heat flux fields for 
the period before and during the pre-ERICA storm are within the North and West 
Atlantic and occurred in the period ranging from 2230 UTC 23 January 1988 to 
1110 UTC 26 January 1988. 

The SSM/I instrument differs from its predecessor (the Scanning Multichannel 
Microwave Radiometer, SMMR) in that the SSM/I does not have a frequency 
from which sea surface temperature can be estimated accurately. Since sea surface 
temperature is required in order to calculate the moisture flux using the bulk 
parameterization method, an independent data source is needed. Two sea surface 
temperature analyses were used to compute the moisture flux fields for the periods 
before and during the pre-ERICA storm. 

The first is a product of the Oceanographic Analysis Branch of the National 
Oceanographic and Atmospheric Administration (NOAA) and is a detailed analysis 
of SST’s in the vicinity of the Gulf Stream Front (GSF). The second is a prod- 
uct consisting of contoured satellite (Advanced Very High Resolution Radiometer, 
AVHRR) SST observations collected over the course of a week. 

The SST field used for the case study is a synthesis of both SST products. The 
primary SST sources were the January 24 and 25, 1988 Oceanographic Analysis 
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GSF analyses. The AVHRR SST map for the week of January 25, 1988 was used 
to fill in the ocean regions not covered by the GSF analyses in in the Northwestern 
Atlantic. 
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3. Methods of Analysis 

The fact that the marine atmospheric boundary layer contains 50% or more of 
the integrated water vapor led Liu (1984) to attempt to estimate monthly mean 
Qa over the North Atlantic ocean using data from the SMMR on the Seasat satel- 
lite. Liu found that the root mean square (RMS) difference between the observed 
monthly Qa and his value computed using a best fit polynomial was 0.78 g kg 1 
and the correlation was 0.99. He makes clear that using time scales shorter than 
a month will result in increased error in Qa using IWV as a predictor. Hsu and 
Blanchard (1989) tested the use of Liu’s fifth-order polynomial relationship be- 
tween monthly IWV and Qa (Liu, 1986) for individual soundings and concluded 
that Liu’s fifth-order polynomial relationship is applicable for individual soundings 
and for soundings located in regions close to coastlines. 

To derive a moisture flux field, the bulk parameterization method is applied in 
the form: 


E = p x Ce x Uio x (Qs - Qa) (1) 

where E is moisture flux [kg (m 2 s) -1 ], p is surface air density, 1.275 [kg m 3 ], Ce 
is water vapor transfer coefficient, U\q is wind speed at a ten meter height [ms 1 J, 
Qs is saturation specific humidity at sea surface temperature [kg kg ’], and Qa is 
specific humidity at a ten meter height [fc<7 

The SSM/I wind speed algorithm is calibrated to give the wind speed at a 
19.5 meter height. We use the Liu, Katsaros, and Businger (1979) [LKB] marine 
atmospheric surface layer model to correct the wind to a 10 meter height. A 
statistical fit is used to derive a (Qs-Qa) value given an SSM/I-IWV value and an 
independent SST value. The LKB model is used to estimate a Ce value for wind 
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speeds less than 10 m s~ l . For wind speeds that exceed 10 m s -1 , we use a Ce 
value of 1.3 x 10 -3 determined by DeCosmo et al. (1990). The p value is assumed 
to be a constant, which results in a maximum density error of ± 0.06 kg m~ 3 . The 
latent heat flux is simply: 


Qe = Lv x E (2) 

where Qe is latent heat flux [W m ~ 2 ], Ly is 4.1868 x (597.31 — 0.56525 x T) x 
(1.0 x 10 3 ) [J fcg -1 ], T is air temperature [C], and E is moisture flux [kg (m 2 s) -1 ]. 

The LKB model takes stratification into account in its iterations. Stratifica- 
tion depends on the air-sea temperature difference and the friction velocity (U,). 
Since no data source exists that gives air temperature at 10 meters (Tio) over the 
spatial coverage of the SSM/I, a stand-in value for Tio must therefore be found. 
A constant temperature lapse rate was assumed where Tio is merely the height 
(10 meters) above the ocean within the surface layer multiplied by the lapse rate 
subtracted from the sea surface temperature. Actual surface layer lapse rates have 
a large variability depending upon origin of the air mass and modification of the air 
mass by the ocean. If the assumed temperature lapse is an overestimate, the corre- 
sponding air temperature at 10 meters will be an underestimate and the resulting 
relative humidity may exceed 100%. 
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4. Evaluation of Statistical Methods 

The best polynomial fits discussed in this study represent the fits that minimize 
the RMS difference. The RMS difference is defined as: 

N 

RMSd=[(l/N)Y,(F l -O i f)i (5) 

1 = 1 

where RMSd is root mean square difference, N is number of independent data 
points, F{ is estimated value, and Oi is observed value. 

Traditionally, IWV has been used to estimate Qa which was then used along 
with a computed Qs (assuming 100% relative humidity at the ocean surface) as 
input for the moisture flux. A more direct method should be one which estimates 
the specific humidity difference (Qs-Qa) between the surface and 10 meter height 
since it determines the direction of the moisture flux. A reasonable step was to 
try a two variable regression of IWV and SST against (Qs-Qa) [hereafter referred 
to as the IWV, SST method]. A polynomial of the form: 

(Qs — Qa) = a x IWV + b x IWV 2 + cx SST 4 -dx SST 2 (6) 

was fit to the January 1988 observations. The units of SST, IWV, and (Qs- 
Qa) are C, kg m -2 , and g kg ~ x , respectively. Coefficients a, 6, c, and d are 
values which minimize the RMS difference. Values for a, 6, c, and d are —0.11676, 
-0.97811 x 10 — 3 , 0.33441, and 0.56958 x 10 -2 , respectively. Although the SST 
(and surface pressure) values are used directly to compute Qs in the development 
of the regression equation, the hope is that SST, in addition to IWV, will account 
for some of the variance observed in Qa and result in a more accurately estimated 
Qs and Qa difference. 
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Figure 1 is a comparison of (Qs-Qa) computed using (6) to the magnitudes of 
(Qs-Qa) from ship data for 124 points. The line shown has a slope equal to 1.0. 
The standard deviation of the observed S/I (Qs-Qa) quantity is 3.73 g kg -1 . The 
standard deviation of the estimated (Qs-Qa) quantity is 3.31 g kg -1 . 

Of the statistical fits tested in Miller (1990), the one employing (6) to estimate 
(Qs-Qa) yielded the smallest RMS differences for both the ‘all’ and ‘no rain’ data 
sets. The difference in the RMS difference between the sample containing ‘all’ 
points and ‘no rain’ points was also smallest for this method. This implies that 
the accuracy of the IWV,SST method is less sensitive to the presence of rain than 
is the IWV versus Qa method (for further details see Miller, 1990). It is apparent 
in Table 1 that the IWV,SST method begins to give useful estimates for (Qs-Qa) 
when the difference exceeds 6.0 g kg -1 . At these magnitudes the RMS difference 
is no more than 23%. This range is also the one where evaporation rates are high. 

The validity of the IWV,SST method could be tested independently given suf- 
ficient information of moisture in the horizontal and vertical dimensions via a 
moisture budget. The evaporation rate could be computed as a residual given 
observed values of moisture storage in a column of the atmosphere, net horizontal 
moisture convergence into this column and an area- averaged precipitation ratefor 
the column. Such a test was not possible at the time of this study since there was 
insufficient wind and moisture data to compute reliable net horizontal moisture 
convergence estimates. 

The ‘no rain 1 points that fell a relatively large distance from the line of perfect 
agreement in Fig. 1 were inspected. These points were in the vicinity of atmo- 
spheric fronts and next to coastlines in cases of offshore flow. Such conditions 
represent cases where a shallow moist layer has formed underneath a dry layer 
of air which is a significant departure from the smooth ‘climatological’ specific 
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FIG. 1. Comparison of (Qs-Qa) based on SSM/I IWV and independent SST analyses 
to Jan. 1988 S/I observations for 124 points with a 2 n< * order fit for (Qs-Qa). 
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humidity profile assumed by applying a statistical method to infer (Qs-Qa) from 
IWV and SST values. Perhaps the most striking outliers are the two with the S/I 
(Qs-Qa) values of approximately -1.0 g kg -1 and 0.0 g kg -1 . These points lie on 
the warm side of a significant horizontal IWV gradient associated with a surface 
cold front. The IWV, SST method is best applied in regions with no large horizon- 
tal IWV gradient. In earlier studies at the University of Washington (McMurdie 
and Katsaros,1985; Katsaros et al., 1989a, b), we have developed a method to objec- 
tively identify frontal zones using IWV maps. IWV gradient values exceeding some 
critical value (e.g. 0.1 (kg m -2 )/ km; Katsaros et al., 1989b) are flagged as being 
located within frontal zones. This method could be used to flag regions within an 
SSM/I overpass where the IWV, SST method should not be applied. However, it 
was not done for this study and, as a result, surface LHF estimation in the vicinity 
of surface weather fronts shown on Figs. 3d, 4d and 8d should be viewed with 
some skepticism. A future study is needed to define the horizontal distance from 
surface fronts where the IWV, SST method becomes applicable. 
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5. The Pre-ERICA Case Study 

A midlatitude cyclone formed off the coast of North Carolina on January 25, 
1988 and deepened rapidly within the following 48 hours. The focus of the syn- 
optic discussion will be on the conditions prior to and during the period of rapid 
deepening. 

a. The Surface Weather Maps 

The surface weather maps used in Figs. 2a and 2b were taken from the National 
Meteorological Center (NMC) Northern Hemisphere Surface Chart analyses. 

At 0000 UTC on January 25 (Fig. 2a), a high pressure cell dominates the 
Northwest Atlantic. At the surface, light winds are observed over the ocean and 
cold temperatures are reported north of the cell center. The cold and stationary 
fronts between 25° and 35° N are attached to a cyclone which has moved east of our 
picture. An open wave cyclone approaches the east coast and has a central pres- 
sure of 1000 mb. By 1200 UTC on January 25 (Fig. 2b), the high pressure cell has 
strengthened to 1031 mb and moved further east. Buoy 44004 (38.5 N,70.6 W) 
reports a surface temperature of 13°C. Chatham, MA reports a surface temper- 
ature of only 4°C so a strong horizontal temperature gradient (0.026°C km -1 ) 
is evident which is a reflection of the influence of the Gulf Stream. The cyclone 
approaching the coast has veered northward and is followed by another cyclone 
located near Chicago, IL. A coastal trough has developed near Cape Hatteras, 
NC. which spawns the eventual rapidly deepening cyclone. Curvature in the 1020 
mb isobar contour at 40°N is indicative of a feature apparently observed by the 
map analyst on satellite imagery. 

The surface analyses used in the following two time periods come from the 
Neiman et al. (1990) article. Much has happened by 0000 UTC on January 26 


15 



Jyj 18 


(a) 0000 irrc 25 January 198* 


(M 1200 UTf 25 January 19*8 


y (c) 0000 UTC 26 January 1988 
* nf !012 1018 


l £ / (d) 1200 UTC 26 January 1988 


F1C. 2. Mean sea level pressure (mb) analyses at (a) 0000 UTC 25 January, (b) 1200 UTC 25 January, (c) 0000 UTC 26 January and (d) 1200 
UTC 26 January 1088. Observed temperatures in (a) and (b) are shown for stations noted in the text. Observed temperatures and pressures in 
(r) and (d) are noted in the upper left and right portions of the station plots, respectively. Analyses in (a) and (b) are taken from NMC Northern 
Hemisphere Surface Charts and the analyses in (c) and (d) are from Neiman et al. (1990). 
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(Fig. 2c). The eventual rapid deepener has formed near Cape Hatteras and moved 
northward as an open wave cyclone. The only other cyclone in the vicinity is 
the one previously located near Chicago. The high pressure cell has stalled to 
the east and weakened. At 1200 UTC on January 26 the storm has continued to 
deepen and moved northward and is on the verge of making landfall (Fig. 2d). It 
experiences a large deepening rate of almost 2 mb hr * from 0000 UTC to 1200 
UTC on January 26 (Neiman et al. 1990). 

b. The SSM/I Derived Fields 

A detailed look at IWV,SST estimated surface latent heat flux fields for some 
SSM/I overpasses occurring both before and during the pre- ERICA case will assist 
in further evaluating the usefulness of our method for obtaining instantaneous flux 
fields over the ocean. 

Several methods and assumptions were utilized for ‘rain-flagged’ data point 
locations outside of frontal zones. Those data points without a meaningful wind 
speed were assigned a wind speed by interpolating values from surrounding SSM/I 
'no rain’ points. Although this interpolation might introduce significant error in 
zones of large horizontal wind shear, it remains as a preferable alternative when 
compared to the probable error introduced by using climatological near surface 
wind speed estimates. 

It was generally assumed that rainrates outside of frontal zones were small 
enough so that a meaningful SSM/I IWV value could be estimated at ‘rain-flagged’ 
locations. Also, the horizontal scale of adjacent ‘rain-flagged data pixels outside 
of frontal zones was generally small enough so that such regions were smoothed 
out in the process of interpolating SSM/I data to the 0.5 grid. Application of 
the IWV,SST method at ‘rain-flagged’ points might be questionable. However, it 
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was concluded previously that the accuracy of the IWV ,SST method for estimating 
(Qs-Qa) is not strongly affected by the presence of rain since the difference between 
the RMS difference for the ‘all’ and for the ‘no rain’ samples was relatively small. 

1 ) Sensitivity to Surface Layer Lapse Rate, Surface Pressure, SST and IWV 
Analyses 

Surface latent heat flux fields were computed for a block of an SSM/I overpass 
to test their sensitivity to large variations in the surface layer lapse rate, surface 
pressure, SST and IWV. In testing the sensitivity of each parameter, surface LHF 
fields were computed with all other parameters held constant as the test parameter 
varied from one extreme to another. There were two trials each for the surface layer 
lapse rate, the surface pressure, SST and IWV sensitivity tests. These tests are 
summarized in Table 2. | A/f | and | AL| represent the largest differences in absolute 
magnitude between corresponding surface LHF maxima and minima, respectively, 
for the various sensitivity trials of each input parameter. The surface LHF fields 
were found to be not sensitive to the assumed surface pressure analysis or the 
surface layer lapse rate for areas where U io exceeds 7ms' but were sensitive 
to the SST, IWV, and U i0 fields and to the surface layer lapse rate for areas of 
low wind speed. The direction of the latent heat fluxes was determined by the 
SST and IWV values, whereas the flux magnitudes were most dependent upon 
L'io- However, the SST and IWV fields do make some contribution to the surface 
LHF magnitudes. Cases where low wind speeds, low(high) SST’s, and high(low) 
IWV values are colocated are ones where the latter two fields have a significant 
impact on the surface LHF magnitudes, but the absolute magnitudes of the fluxes 
are small due to the low U\q values. 
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The apparent insensitivity of the surface LHF patterns to the surface layer 
lapse rate conditions for the pre-ERICA case study makes the use of an assumed 
constant lapse rate reasonable. A surface layer lapse rate of 6.0 °C km -1 will be 
assumed for the remainder of the moisture flux calculations. This value was chosen 
to follow the mean surface to 10 km lapse rate of Nuss and Anthes (1987) and is 
similar to the mean wintertime conditions used by Lau (1978). 

2) A Look at Selected, SSM/I Pre-ERICA Overpasses 

By estimating the fractional error of surface LHF fields derived using the 
IWV,SST method, the usefulness of the fields can be discerned. Since p, Ce , 
t/j 0 , and Ly are identical in surface LHF and surface LHF error calculations, frac- 
tional error is merely a ratio of the RMS difference for the appropriate ‘(Qs-Qa)’ 
bin from Table 1 to the (Qs-Qa) value estimated using (6). The IWV, U\ 9.5, SST, 
surface LHF and surface LHF error calculations for three SSM/I overpasses which 
occurred just prior to and during the existence of the pre-ERICA storm will be 
analysed- 

The IWV maximum east of central Florida for overpass 3091 (Fig. 3a) is to the 
north of the stationary front seen on the 0000 UTC 25 January surface map (Fig. 
2a). The wind speed (Fig. 3b) and local IWV maxima located near Hatteras are 
just east of a strong SST gradient (Fig. 3c) and over water with a temperature near 
22 °C. Clearly the SSM/I has observed a feature not depicted on the 0000 UTC 25 
•January NMC surface map. At the same positions as these maxima are moderate 
upward surface LHF maxima (Fig. 3d), the largest being near Cape Hatteras (273 
IV m ~ 2 ). The area covered in this overpass is located southwest of a high pressure 
cell. Winds are most likely transporting air from the ocean toward land so that 
the atmosphere within view of the SSM/I instrument is maritime in origin. The 
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location of the IWV maxima (Fig. 3a) relative to the GSF is similar to model 
results of Warner et al. (1990). In their study, the Penn State/NCAR mesoscale 
model is initialized with calm winds, horizontally uniform vertical structures of 
temperature and moisture, and a uniform surface pressure field in order to study 
Marine Atmospheric Boundary Layer (MABL) circulations forced by the GSF. 
One result is that low-level moisture convergence increases integrated water vapor 
just to the warm side of the north wall of the Gulf Stream. Another result is the 
appearance of surface latent heat fluxes from the ocean in the 200-300 W m 
range over the warmest water near the GSF. These findings are similar to what we 
observe on the maps of overpass 3091 (Fig. 3). The LHF maxima just off South 
and North Carolina are parallel to the analysed coastal trough shown on the 1200 
UTC 25 January analysis. Without better spatial and temporal data it is difficult 
to conclude if the parallel positioning is related or a coincidence. Three of the 
surface LHF maxima (183, 171, and 273 W m -2 ) are within the 0.3 fractional 
error contour (Fig. 3e). The isolated maximum off Maryland has an estimated 
error greater than 40%. 

An IWV bulge has attained a wavelike shape on the IWV image of SSM/I 
overpass 3097 (Fig. 4a). A stationary front runs roughly east-west on the 1200 
UTC 25 January surface map (Fig. 2b) at the flank of the IWV bulge. This IWV 
feature is colocated with the region of curvature in the 1020 mb isobar contour at 
41°N on the surface analysis, so existence of the isobar curvature appears justified 
by SSM/I imagery. The frontal IWV maxima and IWV gradient have decreased 
from the time of the previous SSM/I overpass (not shown). Wind speed contours 
(Fig. 4b) show little organization at the time of overpass 3097. Speeds in the 8.0 
to 12.0 m s -1 range are found in the 1020 mb isobar curvature region at 40°N 
while along the stationary front they are calm or low (0.0 to 8.0 ms ! ). 
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Maxima within the IWV wave (26.1 and 27.0 kg m ~ 2 ) suggest localized mois- 
ture convergence and marks the active region of the pressure trough. The 37 GHz 
polarization difference map for overpass 3097 (Fig. 5) marks it with a significant 
feature at 40°N, 60°W, where differences are less than 50 K (likely a rain region). 
This ‘rain’ mass is just north of a warm SST eddy (Fig. 4c) and a surface LHF 
maximum (Fig. 4d) at 39°N, 62°W. The proximity of the three features is not a 
coincidence. The 2200 UTC 25 January Sable Island (43.9°N, 60.0°W) sounding 
(Fig. 6) shows distinct upper and lower level moist layers. The upper layer marks 
an approaching warm front, the lower one cannot be the warm front according 
to the 0000 UTC 26 January surface analysis (Fig. 2c) since the surface warm 
front is not yet upon Sable Island and the low level moist layer is positioned next 
to the surface. An SSM/I overpass at 2203 UTC (3104) points to the lower level 
moist layer as being related to the ‘rain’ mass. The 37 GHz polarization difference 
map for that overpass (Fig. 7) shows the advected ‘rain’ mass over Sable Island. 
Although the pre-ERICA storm has matured by 2200 UTC 25 January, the cloud 
mass over Sable Island is too fax north to be associated with it. The warm frontal 
cloud mass accompanying the pre-ERICA storm is just south of 40.0°N. Since the 
‘rain’ mass is in the lower troposphere and it formed on the northwest side of a 
strong high pressure cell, it appears that it received its energy from the ocean via 
latent and/or sensible heating. Surface LHF maxima apparent for an SSM/I over- 
pass preceding 3097 (number 3090) and overpass 3097 confirm that latent heating 
at least was a factor in the ‘rain’ mass development. The 252 W m -2 maximum 
is within the 30% error region. 

The pre-ERICA storm has begun its occlusion and part of the cold front is seen 
south and east of Florida on the 3112 IWV image (Fig. 8a). An IWV maximum 
of 42.7 kg m ~ 2 appears at the cold front. An IWV gradient is apparent behind 


22 



FIG. 4. Fields of (a) integrated water vapor (kg 
m~ 2 ), (b) wind speed at 19.5 meters (m j - 1 ), (c) 
sea surface temperature (°C), (d) surface latent heat 
flux ( W m -2 ) and (e) fractional error derived from 
SSM/I overpass 3097 at 0942 UTC 25 January 1988. 
Sea surface temperatures are contoured every 2°C 
and the fractional error contours are drawn for val- 
ues up to 0.4. The surface front position is from the 
NMC 1200 UTC 25 January 1988 surface pressure 
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FIG. 6. Sounding at Sable Island, N.S. for 2200 UTC 25 January 1988. 
Temperature, solid; dewpoint, dashed; dry adiabats, dot dash; moist adiabat, pecked. 
Winds: full barb, 5 m s' 1 ; pennant 25 m s' 1 . 
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FIG. 7. 37 GHz polarization difference (K) field observed for SSM/I overpas 3104 at 
2203 UTC 25 January 1988. Regions with values less than 50 K are where precipitation is 
most likely to be occurring. 
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the cold front and is associated with strong cold air advection. The 10.0 kg m 
IWV contour reaches to the southern tip of Florida and is evidence of the extent 
of dramatic cold air advection occurring in the southeastern states of the U.S. The 
low IWV values seen on Fig. 8a are consistent with the cold surface temperatures 
reported at 1200 UTC 26 January 1988. Wind speeds at the front range from 
4.0 to 12.0 m s -1 (Fig. 8b). The cold air advection is occurring east of North 
Carolina in regions where winds exceed 12.0 ms *. These winds along with high 
SST’s (Fig. 8c) result in large surface LHF maxima (Fig. 8d) clustered parallel 
to the South and North Carolina coasts. A maximum also is positioned north 
of the Bahamas. The IWV 10.6 kg m ~ 2 maximum well behind the cold front 
marks where strong upward moisture fluxes are modifying the arctic air. As this 
air is advected over the warm ocean, moisture is mixed upward from the ocean 
surface along updraft lines of alternating helical roll vortices which lie horizontally 
within the planetary boundary layer (Hein and Brown, 1988). Cloud streets form 
as a result. The existence of cloud streets was confirmed by looking at the high 
resolution 85 GHz SSM/I liquid water image for overpass 3112 (not shown). A 
look at the 37 GHz Polarization Difference map for overpass 3112 (Fig. 9) shows 
the extent of the developing cloud streets that were seen on the 85 GHz SSM/I 
image. The 60.0 K contour outlines the leading edge of the cloud street region. 
The location of the large upward surface LHF region near the coastline at least is 
confirmed by the fact that cloud streets are present to the east and the surface wind 
is directed so that air moves from the high evaporation rate region to the cloud 
street region. Fractional error (Fig. 8e) for almost the entire overpass is less than 
0.4 with a significant portion less than 0.2. The observed surface LHF maxima 
all lie within a region of less than 30% error. The computed surface LHF fields 
valid at 1200 UTC 26 January 1988 from Neiman et al. (1990) are shown in Fig. 
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FIG. 8. Fields of (a) integrated water vapor {kg 
m -2 ), (b) wind speed at 19.5 meters (m j~ t ), (c) 
sea surface temperature (°C), (d) surface latent heat 
fiux {W rn~ 2 ) and (e) fractional error derived from 
SSM/I overpass 31 12 at 1110 UTC 26 January 1988 . 
Sea surface temperatures are contoured every 2°C 
and the fractional error contours are drawn for val- 
ues up to 0 . 4 . The surface front position is from the 
Neiman et al. ( 1990 ) 1200 UTC 25 January 1988 
surface pressure analysis. 















FIG. 9. 37 GHz polarization difference (K) field observed for SSM/I overpass 31 12 at 
1110 UTC 26 January 1988. The 60 K contour indicates the edge of developing 
convective cloud streets. Regions with values less than 50 K are where precipitation is 
most likely to be occurring. 
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10. Recall that overpass 3112 ocurred at 1110 UTC on the same day. Although 
the spatial overlap between the two analyses is poor, the surface LHF magnitudes 
are comparable for coincident areas in Figs. 10 and 8d as is the tendency for the 
maxima to be oriented in a southwest /northeast direction. 
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6. Summary and Conclusions 

A method was tested which utilized IWV and wind speed from the SSM/I 
along with SST data to estimate surface latent heat flux fields for several SSM/I 
overpasses which occurred during the early and middle development stages of a 
mid-latitude storm. A relationship was tested for estimating (Qs-Qa) as a function 
of SSM/I IWV and an independent SST analysis which indicated some promise 
since it yielded a relatively low RMS difference of 1.96 g kg -1 and was less sensitive 
to the presence of rain than was the other tested method (Qa as a function of IWV). 

Sensitivity of the surface LHF fields to surface layer lapse rate, surface pressure, 
SST and IWV amounts was tested with the conclusion that it was low for surface 
pressure and for the surface layer lapse rate when wind speeds exceeded 7 m s . 
The surface LHF, IWV, SST. U 19 . 5 , and fractional error maps were presented 
for three SSM/I overpasses during the period leading up to rapid deepening of 
the pre- ERICA storm. In eight pre-ERICA SSM/I overpasses that were studied 
(Miller, 1990), results showed that surface LHF fields generally had estimated 
errors exceeding 40% north of 40° N. In these regions, cold air generally is situated 
over cold water so that the saturation vapor pressures at different levels within the 
lowest 10 meters of the atmosphere are nearly equivalent. As a result, evaporation 
is small and the estimated error of the statistical method is large. Most surface 
LHF maxima were colocated with regions of at most 30% error and some regions 
with less than 20% error. These results were computed so that surface wind 
speed error was cancelled out in the calculations. Since it is not well understood 
to what accuracy the surface LHF fields are needed to successfully model rapid 
cvclogenesis, the impact of surface LHF fields estimated from SSM/I and SST data 
on the prediction of rapid cyclogenesis awaits further study. 
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If such fields were to prove useful, employing remotely sensed data for esti- 
mating surface LHF may have most important applications in regions far from 
coastlines and shipping routes. The southern hemisphere oceanic areas are partic- 
ularly data sparse so that methods based on remotely sensed data could broaden 
the tinderstanding of storm development which occurs in those regions. Applica- 
tion of the presented statistical method in regions outside the domain of the case 
study (January in the North Atlantic Ocean) would require an expanded study to 
make the method more general. 

In conclusion, the statistical method developed in this study could provide 
an important contribution to the understanding of rapidly intensifying cyclones 
since the surface latent heat fluxes have been shown (Kuo and Low-Nam,1990; 
Kuo et al. ,1991b) to pre-condition the environment in a way that is favorable for 
cyclogenesis. 
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Table 1: Statistics of ( Qs - Qa)„ t using the second order (Qs - Qa)„ t = f{IWV,SST } polynomial and ( Qs - Qa)„ j, for January 
1088. {NOTE: dQa ok , = (Qs - C?a) e4 ,] 


Bin number: 

i 

2 

3 

Overall RMS 
Difference 
(g/kg) 

Bin interval (g/kg): 

-4.0< dQa ob9 <2.0 

2.0< dQa obs <8.0 

8.0< dQa obs <14.0 

All data points: 

RMS difference {g/kg): 

No. of points = 27 
2.28 

No. of points = 62 
1.80 

No. of points = 35 
2.12 

2.01 

No rain flagged data points: 
RMS difference (g/kg): 

No. of points = 20 
2.26 

No. of points = 49 
1.92 

No. of points = 27 
1.80 

1.96 
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Table 2: Sensitivity trials of four different surface LHF input parameters for a 
mean wind speed greater than 7 m/s. Sensitivity is defined in terms of the largest 
change in maximum (|Ai7|) and minimum (jATj) LHF within the field. 


Parameter 

Trial Values 

Sensitivity 


of Parameter 


Surface Layer 

-6.0 

|AH| = 2 W/m 2 

Lapse Rate 

9.8 

|AL| = 1 W/m 2 

(C km ~ ] ) 



Surface Pressure 

940 

|Ai/| = 1 W/m 2 

(mb) 

1050 

\AL\ = 0 W/m 2 

SST 

0.0 

\AH\ = 900 W/m 2 

(C) 

27.0 

|AL| = 500 W/m 2 

IWV 

0.0 

|Ai?| = 550 W/m 2 

(kg m~ 2 ) 

50.0 

\AL\ = 450 W/m 2 
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ABSTRACT 

Until recently, the scarcity of meteorological observations over polar areas has limited 
studies of high-latitude weather systems, but now data from polar orbiting satellites offer a 
new opportunity to observe and describe these systems. TOVS data have been used 
successfully for delineating synoptic and sub-synoptic systems since they provide the 
vertical temperature structure of the atmosphere: SSM/I observations have proved valuable 
for analysing storms through water vapor and rain determinations. These positive results 
prompted us to analyze simultaneous TOVS and'SSM/I observations obtained during a cold 
air outbreak over the Norwegian Sea. After a description of the instruments and the 
retrieval schemes, the mutually supporting information from these two independent 
instruments is discussed. Implications for the monitoring of polar lows are presented. 
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1. Introduction 


Until recently, the scarcity of meteorological observations in polar regions has seriously 
hindered the study of weather systems in these areas. The advent of polar orbiting satellites 
with the capability of retrieving quantitative meteorological information now offers a new 
opportunity. Because of the large overlap at high-latitudes of consecutive swaths of data 
from polar orbiting satellites, one can obtain frequent samples of mesoscale weather 
systems and follow their development and decay. 

From the TIROS-N Operational Vertical Sounder (TOVS) aboard NOAA satellites, it is 
possible to extract the vertical structure of the atmosphere in terms of temperature, moisture 
and cloudiness. Atmospheric temperature profiles can be obtained with good accuracy in 
most weather conditions (Le Marshall, 1988). Unfortunately, the determination of other 
parameters such as water vapor profiles or total water vapor content is impossible if the 
cloudiness is greater than a certain threshold; therefore, in polar regions where cloudiness 
is prevalent, such determinations often fail. 

One way of alleviating this problem over the sea is to use the observations of the Special 
Sensor Microwave /Imager (SSM/I) launched in June 1987 on a satellite in the United 
States Defense Meteorological Satellite Program (DMSP) (Hollinger et al., 1987). Several 
studies have already shown that integrated atmospheric water vapor (precipitable water), as 
well as surface wind speeds can be retrieved over the open ocean from this instrument with 
good accuracy (Alishouse et al., 1990 ; Petty, 1990 ; Goodberlet et al., 1989). 

In this paper, we present the application of both TOVS and SSM/I observations to the 
analysis of a cold air outbreak over the Norwegian Sea between 29 February and 1 March 
1988. After some background on both instruments and the retrieval algorithms in section 2, 
section 3 follows with an overview of the synoptic situation. The TOVS and SSM/I 
retrieved parameters for the case of interest are then presented in sections 4 and 5 
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respectively. Finally, in section 6, the complementary nature of the information from these 
two types of satellite data is discussed. We show that the combination can be of 
considerable value for the monitoring of polar lows, which often occur in these regions 
during wintertime and which are difficult to forecast due to their small size and short 
lifetime. 

2. The sensors and retrieval algorithms 

2. 1 . The TOVS instrument 

2.1.1. Description. The TIROS-N Operational Vertical Sounder, which flies aboard 
NOAA satellites consists of three passive vertical sounding instruments (Smith et 
al.. 1979): the High resolution Infrared Radiation Sounder (HIRS-2), a radiometer with 19 
channels in the infrared and one in the visible: the Microwave Sounding Unit (MSU), a 
microwave radiometer with 4 channels in the vicinity of 55 GHz; and the Stratospheric 
Sounding Unit (SSU), a pressure-modulated infrared radiometer with 3 channels near 15 
urn. 

Table 1 gives the characteristics and primary use of observations made in the various 
spectral channels of the TOVS. The HIRS-2 ones are distributed along the whole infrared 
spectrum (from 3.7 to 15 (im). The ability of the 50 GHZ channels to penetrate clouds 
make them indispensable, since clouds are generally opaque to the infrared channels. 

The spatial resolution of the infrared sounder is about 17 km at nadir whereas the spatial 
resolution of the microwave sounder is 100 km at nadir, the swath width being about 2400 
km. Besides the TOVS instrument, there is also an imager aboard the NOAA satellites, the 
Advanced Very High Resolution Radiometer (AVHRR) with a spatial resolution of 1.1km . 
Normally two TIROS-N satellites operate simultaneously ; at low latitudes, each satellite 
observes a point on the earth every 1 2 hours, but at high latitudes, overlap between swaths 
allows more frequent sampling. 
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2.1.2. The retrieval algorithm. The TOVS retrieval algorithm used in this study is the 
Improved Initialization Inversion (31) algorithm (Chedin and Scott, 1985 and Chedin et al., 
1985), developed at Laboratoire de Meteorologie Dynamique. This method determines the 
following physical parameters : air mass type, temperature profiles (or geopotential 
thicknesses), thermal winds, microwave surface emissivity, cloud heights and amounts, 
total water content and relative humidities in 3 layers. Only HIRS-2 and MSU observations 
are currently used, so that the retrievals are obtained from the surface to a pressure level of 
10 hPa (by adding SSU observations, it would be possible to reach 1 hPa). 

The inversion relies upon a priori knowledge of the parameters (the atmospheric structure) 
and the associated calculated brightness temperatures for each channel of the vertical 
sounder. This a priori knowledge is in the TOVS Initial Guess Retrieval (TTGR) data set, 
created once and for all. Provided the observed brightness temperatures correspond to clear 
areas, or have been properly "cleared" (decontaminated from clouds) by the so-called psi- 
method, the 31 procedure consists of two principal steps: first a search is made for the best 
initial guess solution in the TTGR data set ; finally the solution is obtained by a maximum a 
posteriori probability estimation procedure using the Jacobians associated with the initial 
guess, in the TIGR data set. 

Retrievals of water vapor quantities (relative humidity, integrated water vapor) are obtained 
for 3 layers delimited by the levels 1000, 850, 500 and 300 hPa corresponding 
approximately to the center of HIRS-2 channels 10, 1 1, 12 weighting functions ( Fig. 1 ). 
Following the temperature inversion, predicted brightness temperatures associated with the 
initial guess are corrected for the deviations between the initial and final temperature 
profiles. These new calculated brightness temperatures are then used as initial guess for 
water vapor and surface temperature retrievals. The method is based upon a simultaneous 
physical inversion of water vapor amounts and surface temperature ( HIRS-2 channel 10 is 
sensitive to the surface temperature, and use is also made of HIRS-2 window-channels 8 
and 18 ). Reliable water vapor retrievals become impossible when cloud coverage is large, 
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Fig. 1 : Weighting functions of the 6.3 |im ("water vapor ’) and 1 1 M- m ( windows ) 
HIRS-2 channels ( from Smith et al., 1979 ). 
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i.e. greater than 60 %. In such cases, holes due to rejection by the algorithm appear in the 
retrieved water vapor fields (see e.g. Tahani et al„ 1989). 

All parameters are retrieved at a spatial resolution of 100 km. Quality control is 
incorporated into the procedure in order to detect and to reject either bad satellite data or bad 
retrievals. For this purpose, all retrievals pass a series of tests using the deviations 
between observed, initialized and retrieved profiles (Chedin, 1988). In cases with 
precipitation detected by MSU observations, retrievals must also be rejected (Chedin et al., 
1989). In the arctic and subarctic environment, temperature profiles are obtained with good 
accuracy and most of the time compare well to data from radiosondes (Claud et al., 1991). 

2.2. The SSMil instrument 

2.2.1. Description. The Special Sensor Microwave/Imager (SSM/I) was launched by the 
Defense Meteorological Satellite Program (DMSP) on June 19, 1987. This seven channel 
microwave radiometer measures upwelling polarized microwave radiation from the 
earth/atmosphere system at four frequencies : 19.35, 22.2, 37, and 85.5 GHz. Two 
polarizations, vertical (V) and horizontal (H), are available for each frequency, except for 
22.2 GHz which has only vertical polarization. The spatial resolution ranges from 55 km at 
19 GHz down to 15 km at 85 GHz. Further details concerning the instrument are given by 
Hollinger et al., (1987). 

Radiation received by the SSM/I comes from three sources : atmospheric emission in the 
direction of the satellite, surface emission, and downward atmospheric emission reflected 
from the ocean surface. Changes in the total microwave radiance observed by a given 
channel are dominated by variations in the atmospheric transmittance and, to a lesser extent, 
by changes in the ocean surface emissivity due to wind-induced roughening and foam 
production. At SSM/I frequencies, the most important atmospheric attenuators/emittors are 
oxygen, water vapor and liquid water (Fig. 2). 
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Fig. 2 : Illustration of the positions of SSM/I channels relative to oxygen and water vapor 
lines in the microwave spectrum ( 0 - 240 GHz), after Waters, 1976. The zenith opacity 
refers to the integrated column of a standard atmosphere. 


48 



Because each of the seven channels of the SSM/I responds differently to changes in 
integrated water vapor, integrated liquid water, surface wind speed, and certain "minor" 
parameters, it is possible to combine channel brightness temperatures in a way that is 
sensitive to a single parameter of interest and approximately orthogonal to the rest. The 
general problem of optimally retrieving environmental parameters from the SSM/I has 
been discussed in some detail by Petty ( 1990). 

2.2.2. The retrieval algorithms. Here we confine our attention to two parameters which 
are known to be measurable by the SSM/I with reasonable reliability : the integrated water 
vapor content and the surface wind speed. These are both obtained with a spatial resolution 
of about 50 km, though the sampling density is somewhat higher. Because of the relatively 
small amounts of water vapor and liquid water encountered in polar air masses, simple 
linear water vapor and wind speed algorithms are quite satisfactory for the purposes of this 
study. 

We have developed an algorithm for integrated water vapor specifically intended for use in 
the polar oceanic regions considered in this study; it depends simply on the difference 
between the 22V and 19V channel brightness temperatures : 

IWV (kg/m 2 ) = 0.98 (TB22V-TB19V) -3.5 (1) 

where TB19V stands for the 19 GHz brightness temperature in vertical polarization, and 
TB22V, 22 GHz brightness temperature in vertical polarization. 

The coefficients for the above relationship were derived from the analytic brightness 
temperature model of Petty ( 1 990), assuming atmospheric water vapor and temperature 
profiles characteristic of air masses with less than 20 kg/m 2 of total water vapor, as 
determined from the global radiosonde statistics presented by Petty (1990). An additional 
assumption was that the mean surface wind speed was near 10 m/s and that cloud water in 


the amount of approximately 0.1 kg/m 2 was present and had a temperature near 270 K. 
Since the brightness temperature difference appearing in (1) is only weakly sensitive to any 
of these assumptions, we expect that maximum uncertainties in retrieved IWV do not 
exceed 1-2 kg/m 2 for the cases considered and are probably less, provided only that IWV 
does not greatly exceed 15 kg/m 2 . For higher values of IWV, the assumption of linearity 
fails, and (1) will give an increasingly serious underestimate of the true water vapor 
amount. 

For surface wind speed determination, we used an adaptation of the algorithm of 
Goodberlet et al„ (1989), who derived the following coefficients via a "weighted linear 
regression" of 708 buoy wind speeds against the coincident SSM/I brightness temperatures 
from the 19V. 22V, 37V, and 37H channels: 

SW (m/s)= 147.90+1.0969 TB19V-0.4555 TB22V-1.76 TB37V+0.786 TB37H (2) 

SW is the wind speed referenced to a height of 19.5 m above the surface. Our own 
extension of this algorithm includes spatial smoothing and slightly modified criteria for 
flagging regions contaminated by rainfall. These minor modifications are described more 
fully by Petty (1990). 

Under rain free conditions, a standard error of about 2 m/s was reported by Goodberlet et 
al., ( 1989) for winds in the range 3-25 m/s. This error statistics may not be representative 
for the regions and extreme meteorological conditions considered in the present study. 
However, our own experience suggests that the retrieved values are reasonable, though 
with an apparent tendency to underestimate wind speeds above 25 m/s. 

For completeness, we also examined two other parameters which can be inferred from 
SSM/I data, the cloud liquid water content and an index of scattering by large ice particles. 
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This index is indicative of convection. Large ice particles, graupei or hail, are required for a 
substantial signal. The algorithms used have been validated more qualitatively than 
quantitatively. For a complete description, the reader is again referred to Petty (1990). 

3. Synoptic description of the 29 February to 1 March 1988 situation 

Maps of the European Meteorological Bulletin (EMB, see references) as well as AVHRR 
images have been used for describing the synoptic situation. The three day period from 29 
February to 1 March 1988 was characterized by a cold air outbreak which had abruptly 
begun over the Norwegian and North Seas on February 28. The northerly flow associated 
with the cold air outbreak was reinforced during this period by a strong east-west pressure 
gradient accompanying a cyclone over Denmark and southern Sweden with a central 
pressure 980 hPa and a strong anticyclone situated south of Greenland with a central 
pressure of 1050 hPa on 28 February at 1200 UTC (Fig. 3a ) . An upper level disturbance 
(500 hPa) was associated with the cold air region delimited by the fronts at the surface. 

During February 29, the region covered by cold air expanded southward (Fig. 3b ); the 
front extending from northern Norway to Novaya Zemlva evolved into a low center in the 
Barents Sea on the 1st of March. A stationary front lay west of Great Britain. Surface 
winds were quite strong in frontal areas as well as in the North Sea. At 1200 UTC, almost 
the entire area was cloudy except a region close to Greenland. The 500 hPa disturbance 
moved southward together with the cold air at the surface. 

Fig. 4 displays the AVHRR image in channel 4 at about 0900 UTC ; it indicates the 
presence of a frontal cirrus band between Iceland and the British Isles. Another cirrus band 
is visible over Scandinavia. Most of the region corresponding to the strong pressure 
aradient is covered with cumulonimbus and well-defined cloud streets can be seen west of 
Spitzbergen. Note also the formation of a vortex at about 76 °N, 13 °E. 
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On March 1 at 1200 UTC ( Fig. 3c ), a col developed northeast of Iceland in a trough 
connecting two lows, one over Iceland, the other close to Bear Island. The cold front, 
which by now had reached the African coast, was weakening. The area of warm air over 
the Atlantic Ocean had moved eastward and the winds over the Norwegian and North Seas 
were much weaker than on the previous day. 

The AVHRR image in channel 4 at 0845 UTC ( not shown ) gave further indications of the 
evolution over the previous 24 hr : both the cirrus band over Scandinavia and the one 
between the British Isles and Iceland were seen to have moved eastward, and while the 
cloud streets west of Spitzbergen were less dramatic, the small vortex first identified on the 
29th had become much more pronounced and was bounded on the east by a cloud-free 
area. It is likely that this system corresponds to a polar low identified at the Norwegian 
Meteorological Institute (DNMD at 1200 UTC (position : 74 °N, 18 °E). 

4. Results of the 31 algorithm 

Three NO A A- 10 passes covering the area were obtained from the Norwegian 
Meteorological Institute ( DNMI ) for 29 February at 0856 UTC, 29 February at 1707 
UTC. and l March at 0833 UTC. Each pass is of about 10 minutes duration. Cloud related 
results, geopotential thicknesses and water vapor retrievals are discussed below. 

4.1 . Cloud related results 

Cloud detection results for the 29 February pass are presented in Fig. 5 (these correspond 
to the AVHRR image in Fig. 4). The situation is found to be quite cloudy by the 31 
algorithm, as would be expected based on the AVHRR image. The high cloud systems 
over Iceland and Scandinavia have been well detected. The 31 cloud top pressure is in the 
range 300-400 hPa, which seems reasonable, though it cannot be independently verified. 
Rather low clouds over the Atlantic Ocean west of the polar front have also been well 
detected, and the top pressure deduced from 31 is in the range 600-800 hPa, which is 
realistic in view of the synoptic situation and the infrared image (Fig. 4). Nevertheless, 
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some clear areas are identified by the 31 system in regions where there are obviously 
cumulonimbus created by surface heating. The failure to detect these clouds is probably due 
to the fact that they do not cover the whole box (cloud detection tests are performed within 
each HIRS-2 spot, but 2 clear spots are sufficient in a box to declare it clear). For the 
following passes, these remarks also apply, the large cloud systems being well detected, 
except again most of the cumulonimbus clouds. 

4.2. Geopotenrial thicknesses 

Fig. 6a-c display the 1000-500 hPa geopotential thicknesses. On February 29 ( Fig. 6a ), 
two general regions of large thickness gradient can be identified : one corresponds to the 
frontal zone extending from Greenland across Iceland to just west of Ireland (Fig. 2a, 
weather map); the other includes the region poleward of 65 °N. The latter region is 
characterized by two north-south axes of lower thickness near 5 °W and 20 °E respectively, 
these features persist throughout the following 24 hr ( Fig. 6 b, c ), and are seen to drift 
eastward by only 5 degrees longitude during this period. 

Unfortunately, the full extent of the cold air cannot be seen in the evening pass on the 
29th, since this pass only covers a portion of the North Sea. However, on the 1st of 
March, the southward advancement of the cold air has been halted : the front between 
Greenland and the British Isles has moved eastward of its position on the previous day and 
the front in the Barents Sea is sharper. 

A striking feature, however, is the large number of rejections, especially for 29 February 
and 1 March passes in the Atlantic Ocean and in the North Sea. In the case of the North 
Sea, these rejections are due to precipitation. Over the Atlantic Ocean, the situations found 
in the TIGR data set were not close enough to the observed situations and were rejected by 
the quality control tests (see section 2.1.2). 

The 31 inferred 1000-500 hPa geopotential thicknesses for these three passes were found to 
generally agree with the analyses (not shown) provided by the ECMWF (European Center 
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a. 29 February 1988, 0856 UTC 

b. 29 February 1988, 1707 UTC 

c. 1 March 1988, 0833 UTC 

Areas of rejections ( section 4.2 ) are shaded grey. The dashed line represents the swath 
edge. 






for Medium range Weather Forecasts) ; similar large-scale structures appear though the cold 
tongues on both sides of Spitzbergen are less pronounced on the ECMWF products. 

The 1000-850 hPa geopotential thickness field (Fig. 7) on the morning of 29 February is 
qualitatively very similar to the corresponding 1000-500 hPa thickness field. Note however 
the presence of a tongue of warm air extending from England to about 66 °N just off the 
western coast of Norway. 

The 31 derived 500-300 hPa geopotential thicknesses displayed in Fig. 8 for the same pass 
;gain show the gradient associated with the front extending from Greenland to west of 
Ireland. 

4. 3. Water vapor 

Total retrieved water vapor content is displayed in Fig. 9a-b for the morning of 29 
February and 1 March, respectively : in several regions, no determination has been 
performed, due to either rejections in the temperature and/or water vapor profiles 
calculation or cloudiness larger than 60%. 

For both passes, two regions can be distinguished : the warm air side of the front is seen to 
be generally characterized by water vapor values of 8-12 kg/m 2 , while the cold air mass 
north and east of the British Isles is considerably drier. In this latter region, the warmest air 
in the vicinity of 60 °N, 10 °W is characterized by water vapor content of about 7 kg/m 2 , 
while the coldest air near the ice edge presents water vapor content as low as 1 kg/m 2 . Note 
that the fields of retrieved integrated water vapor values are found to be qualitatively very 
similar to the 1000-500 hPa geopotential thickness fields, including even the two thickness 
tongues both sides of Spitzbergen described earlier. 

It is also possible to derive the water vapor content for 3 separate layers : 1000-850 ; 850- 
500 ; 500-300 hPa. Most of the water vapor (from 30 to 50 %) is found in the layer 1000- 
850 hPa (see Fig. 10 a-b), with values ranging from 0.5 to 5.5 kg/m 2 . The warm air 
advection west of Norway observed in the thicknesses is also apparent in the water vapor 
content which has a relative maximum of 2.5 kg/m 2 . In the 800-500 hPa layer (Fig. 1 la- 


58 



Fig. 7 : 1000-850 hPa 31 geopotential thicknesses (in decameter), 29 February 1988, 0856 
UTC. Areas of rejections are shaded grey. 
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Fig. 9 : 31 total precipitable water (in kg/m 2 ) 

a. 29 February 1988, 0856 UTC 

b. 1 March 1988, 0833 UTC 
Areas of rejections are shaded grey. 
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b), values range from 1 to 6 kg/m 2 , which represents about 40 to 50 % of the total content, 
but above 500 hPa, values are extremely low ( less than 0.5 kg/m 2 ). 

In relatively moist areas, both the lower and the intermediate layers contain 40 to 50 % of 
the water vapor, whereas a smaller percentage (30 %) is found in the lowest layer in the 
very dry areas. 

If we consider now the evolution from the 1 st orbit to the 3rd one, layer by layer, it appears 
that the temporal evolution in the 1000-850 hPa layer is very weak for almost the whole 
pass, except over the Atlantic Ocean north of Spain, where the values increase from 3.6 to 
5.4 kg/m 2 . The same tendencies appear for the 850-500 hPa layer. 

5. Analysis of SSM/I parameters 

For this period, there were several SSM/I passes covering about the same area as the TOVS 
observations (Table 2). Poleward of 57 °N, successive SSM/I swaths overlap, so that there 
is complete geographic coverage twice daily. Five SSM/I composites, each consisting of 3 
or 4 overlapping swaths (28 February evening, 29 February, early morning and evening, 1 
March, early morning and evening) are considered for presenting the SSM/I results. 

5.1. Water vapor. 

Selected composite SSM/I water retrievals are plotted in Fig. 12, using the same contour 
intervals and map projections as for the TOVS derived products in Fig. 9. In contrast to the 
TOVS, the SSM/I is limited to retrieving integrated water vapor over ice free water ; a 
dashed line is used to indicate the edge of the pack ice beyond which no SSM/I data can be 
used. 

The SSM/T integrated water vapor image on the 28th (not shown) exhibits a characteristic 
water vapor maximum ahead of the polar front, accompanied by a sharp decrease across it. 
Quickly, the dry air advances southward, eastward and westward during the early morning 
of 29 February. The British Isles are at that time entirely in the cold air which has even 
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reached France. Southeast of Iceland, the water vapor content has increased dramatically in 
the warm sector of the incipient wave on the polar front (63 °N, 13 W). 

The SSM/I retrievals for the evening of 29 February (Fig. 12a) and the early morning of 1 
March (Fig. 12b) clearly confirm the presence of a strong dry air incursion over the 
Norwegian Sea, with extremely low water vapor values covering most of the open ocean 
north of 60 °N. Throughout the region covered by the cold air mass, SSM/I retrieved water 
vapor values show a week north-south gradient in water vapor amount, ranging from 1-2 
kg/m 2 over the open water just south of Spitzbergen up to about 6-7 kg/m 2 near the 
advancing edge of the front extending from Greenland to the British Isles. 

The SSM/I retrievals also clearly reveal the sharp boundary between the arctic air and the 
somewhat moister air mass to the south and west of Iceland and the British Isles, 
respectively. Based on the work of Katsaros et al. (1989) and McMurdie and Katsaros 
(1985), who found that a surface cold front is typically located on the moist edge of the 
strongest gradient in the integrated water vapor field, these maps suggest a surface frontal 
position extending from southern Iceland to nonhem Ireland at around 2000 UTC 29 
February (Fig. 12a) and then from the northeastern tip of Iceland to central Ireland near 
0530 UTC 1 March (Fig. 12b). These positions are in fact found to be generally consistent 
with the operational synoptic analyses presented for 1200 UTC on 29 February and 1 
March ( Fig. 3 b and 3 c). 

On the warm side of the front, the SSM/I data also show a relatively narrow water vapor 
maximum of 12-14 kg/m 2 ; such a maximum is commonly observed in this position near a 
surface front and is believed to result from a combination of enhanced convergence and 
potential advection of moist air along the frontal zone. 

Throughout this period, the Barents Sea (characterized by a low surface pressure), remains 
relatively moist (4 to 6 kg/m 2 ), compared to the region of the cold air outbreak. 
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5.2. Winds. 


On the 28th of February, SSM/I surface winds (not shown.) are strong in the cold air region 
(greater than 16 m/s), but in a large area over the North and Norwegian Seas, they could 
not be determined because of precipitation. There is a sharp gradient in the winds across the 
two frontal zones, one extending from southeast of Iceland to the British Isles, the other 
from northern Norway to Novaya Zemlya. These strong winds correspond to strong 
pressure gradients appearing in the EMB map at 1200 UTC, where wind speeds of 30 m/s 
are reported in the southern part of the Norwegian Sea. However, some SSM/I winds in 
the Greenland Sea are higher than suggested by the EMB pressure fields. A similar feature 
has been observed with a different algorithm in mid-latitude and can be explained by the 
fact that SSM/I wind speeds would be sensitive to surface roughness variations due for 
example to atmospheric stratification or sea state not directly related to th. * wind (swell, 
current effects). 1 

Values in the range 4-12 m/s occur in the westernmost moist areas and seem in fair 
agreement with EMB wind repons at 1200 UTC. 

The situation during the early morning of 29 February (not shown) does not show much 
evolution since the previous day, with still strong gradients across the frontal zones. 
However, less precipitation is found in the Norwegian Sea, but more in the North Sea. An 
. area near 65 C N, 10 °W is found to have lower wind speeds (4-8 m/s) than some hours 
earlier, when they were larger than 14 m/s. These values are difficult to validate, since there 
are no surface reports in this area at that time (0555 UTC). However, at 1200 UTC, EMB 
analysis reports low wind speeds and slack pressure gradients. 

On the 29th of February, evening overpasses (Fig. 13a), precipitation occurs only in the 
frontal area between Iceland and the British Isles, so the winds can be determined for 
almost the whole swaths. The contrast between wind speeds in the cold air and in the warm 
air tend to be weaker. Largest SSM/I wind speed values (in the range 10-14 m/s) are still 
found in the area of cold air (Norwegian and North Seas), but also west of the front 
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Fig. 13 : same as Fig. 12, but for surface wind speeds, derived from Eq..2 and only for 
evening overpasses of 29 February (a) and 1 March (b). Values range from 0 to 20 m/s. 
The absence of determination ( grey areas over sea ) indicates the presence of precipitation. 
The dashed line represents the approximate sea ice edge. 
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extending from Iceland to the British Isles. SSM/I wind speed values in the Barents Sea ( 8 
to 12 m/s) are greater than what one might infer from the pressure patterns appearing in the 
EMB maps (Fig. 3b, 29 February, 1200 UTC). Reasonable agreement is, however, found 
west of the Barents Sea between SSM/I wind speeds and wind repons ( e.g. 8 m/s at 
75 °N, 20 °E. and 10 m/s at 76 °N and 25 °E). These differences in agreement may be due 
to differences in stratification ( in the Barents Sea, the southerly winds bring warm air over 
the ocean, while west of this area, cold air streams off the arctic sea ice over warmer 
water), which results in different efficiency of momentum transfer through the atmospheric 
boundary layer. 

On the 1st of March early morning, SSM/I wind speeds (not shown) in the Barents Sea are 
in better agreement with the EMB analysis.The area with lower wind speeds southeast of 
Iceland has extended eastward whereas the area of high wind speeds over the Norwegian 
Sea is smaller. This wind distribution agrees well with the EMB analysis at 0000 UTC (not 
shown). 

In the evening of the 1st of March (Fig. 13b), rather large values (greater than 12m/s) are 
only found in small zones between Norway and the British Isles ( EMB repons at 1200 
UTC a wind speed of 10 m/s at 59 °N, 25 °W ). 

5.3. Signature of a polar low. 

We have also looked at two other parameters derived from SSM/I, the cloud liquid water 
content and an index of scattering by large ice particles at 85 GHz. The valuable 
complementary information provided by these two parameters will be illustrated with the 
detection of a polar low on March 1, 0401 UTC. In Fig. 14 a which displays cloud liquid 
water at that time, a large region of liquid water coincides with the region of cumulus 
advection already seen in the AVHRR image on 29 February at about 0900 UTC (Fig. 4) 
and a linear feature is identifiable between 72 and 75 °N and at about 17 °E. In Fig. 14 b, 
the scattering index delineates this feature more clearly along with the more intense 
convective clouds close to the Norwegian coast. The detection of this polar low is 
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corroborated by the A VHRR image on the same day at 0845 UTC (not shown) and by the 
DNMI analysis, which identified this low at 1200 UTC at a position of 74 °N, 18 °E. In the 
SSM/I evening passes of March 1 , there is no longer a signature (nor does DNMI list it any 
longer at 1800 UTC), one can surmise that it has either disappeared or is very weak. 

6. Discussion 

We have used both TOVS and SSM/I observations to study a cold air outbreak over the 
Norwegian Sea. Each instrument gives a special view of the weather system: the TOVS 
instrument provides a vertical description of the atmosphere in terms of temperature, 
moisture and cloud cover at a coarse spatial resolution (100 km), whereas the SSM/I 
instrument yields the integrated water vapor content and the surface wind speed at a higher 
spatial resolution (50 km). Moreover, the sensitivity of SSM/I to cloud liquid water and 
scattering by large ice particles makes it a powerful instrument for detecting embedded 
convection. 

Comparing water vapor fields obtained from SSM/I and TOVS, one sees reasonable 
qualitative agreement for large-scale structures, in spite of time difference and of a 
significant number of rejections for TOVS derived water vapor. In regions where the 
paucity of conventional observations makes validation difficult, it is reassuring to have two 
independent instruments produce similar large-scale features. However, SSM/I shows 
much sharper gradients and generally more detailed subsynoptic structure, due largely to its 
better spatial resolution and its more complete coverage ( SSM/I determines integrated 
water vapor content through clouds and rain). Throughout the region covered by the cold 
air mass, both SSM/I and TOVS show a weak north-south gradient in water vapor amount, 
but unlike the TOVS retrievals, the SSM/I ones clearly reveal the sharp gradient between 
the arctic air and the moister air south and west of the polar front. The TOVS retrievals yield 
only marginal evidence for the front related features in the water vapor fields, because of its 



less favorable spatial resolution and the problem of cloud contamination throughout most of 
the vicimcv of the front. 

The 31 water vapor fields are found to be qualitatively very similar to the geopotential 
thickness fields of the same layer ( including even the two tongues described in section 
-i.2). This observation is consistent with the expectation that cold air will generally be drier 
than warm air. On the other hand, the possibility exists that some of the similarities 
between the TOVS thickness and water vapor fields are an artifact of the retrieval. The fact 
that SSM/I water vapor retrievals show significantly less structural similarity to the TOVS 
thickness field would indicate that the second explanation is valid. However, the large 
number of rejections in the TOVS water vapor field may also have contributed to these 
similarities. 

Quantitative comparisons between TOVS and SSM/I integrated water vapor content give 
on average a root mean square difference of 20 to 30 % of the mean TOVS values. Given 
that only a few conventional measurements are obtained over the area of interest, and that 
these measurements may not be representative for large areas, SSM/I and TOVS humidity 
estimates which are of comparable quality to radiosonde measurements could be important 
for meteorological purposes as underlined bv several authors for TOVS (Illari 1989, 
Zhang 1989). These authors mention that the impact of satellite derived humidity 
parameters on humidity analysis are particularly important for data-sparse areas over the 
ocean. 

SSM/I has the advantage of providing integrated water vapor content over open ocean in 
spite of cloudiness and rain (if not too intense). From TOVS data, no water vapor 
determination is performed in presence of rain or if the cloud cover is substantial, but 31 
provides water vapor content over sea ice and land, where SSM/I can not discriminate the 
water vapor signal from the background. There is currently little validation of this 
parameter over sea ice. but the very iow amounts found look reasonable. Over land, 
comparisons between TOVS derived integrated water content and radiosonde or microwave 


ground-based radiometers inferred content show that the difference in the estimation is 
about 20 % rms (Tahani, 1991). The coupling of both TOVS and SSM/I allows one to 
obtain "quasi-continuous" fields of water vapor in space, which are required for numerous 
meteorological applications. 

Moreover, some information about the vertical distribution (3 separate layers : 1000-850, 
850-500, and 500-300 hPa) can be obtained with TOVS. Even coarse vertical resolution 
gives insight into the physical processes involved. 

Surface wind speeds derived from SSM/I, although sometimes stronger than what would 
be expected from conventional reports, appear to be realistic when compared to weather 
maps. It is nevertheless true that regions of a swath contaminated by rain (i.e. most of the 
time, the interesting parts of weather systems) are eliminated, which represents a 
shortcoming. However, regions of high wind speeds far enough from fronts are well 
delineated as shown in this study. 

A polar low has been identified by two different parameters deduced from SSM/I (cloud 
liquid water and scattering by large ice particles). This polar low was undetected by TOVS 
observations and there is no signature of it on the surface analyses because of its small 
spatial scale and its very short lifetime. Further work dedicated to polar lows detection 
confirms this potential of SSM/I observations (to be published). An earlier study mentioned 
the utility of multispectral microwave radiances for analyzing these disturbances (Gloersen 
et al„ 1989), and more particularly the 18 and 37 GHZ frequencies (two frequencies 
avalaible on the Scanning Multichannel Microwave Radiometer onboard the Seasat and 
Nimbus-7 satellites). There is no doubt that the high spatial resolution as well as the unique 
properties of the SSM/I 85 GHz channels added to the 19 and 37 GHz channels makes 
SSM/I a remarkable instrument for further studying polar lows. Coupled with the TOVS 
instrument, the SSM/I increases our capabilities for monitoring polar low developments in 
data-sparse regions. 
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Figure captions 


Table / : Characteristics ofTOVS sounding channels. 


r abk 2 : Ctaracertsttcs of SSMi, souths (date, orbit number and time CTO and 

■corresponding significant meteorological features. 

= ig . l : Weighting functions of the 6.3 pm ("water vapor") and 11 (un ("wmdows ) 
HLIRS-2 channels ( from Smith et al.. 1979 V 

ccvi/T channels relative to oxygen and water vapor 
r . . nincrT-irion of the positions ot SSM/i cnanneis icm j . , . 

F,g. 2 Illustnmon 01 P Waters. 1976 . The zenith opacity 

lines in the microwave spectrum t 0 - -40 OHz), alter 

refers to the integrated column of the 1962 US Standard Atmosphere. 

Fi , 3 : surface chart (tedrawn from the Emopean Meteorological Bulletin) a, .200 UTC 

a. 28 February 1988 

b. 29 February 1988 

c. 1 March 1988 

Hg . a : AVH RR .mage, channel 4. 29 February 1988. 0856 UTC C Courtesy of the Royal 
Dutch Meteorological Institute , [KNMI1). 

R . - ; 31 cloud detection results, 29 February 1988, 0856 UTC. Clear 

bv'dark dots, cloudy ones are shaded grey. Blank zones correspond to cal.bmt.on . 

Fis. 6 : 1000-500 hPa 31 geopotenual thicknesses (in decameter) . 

a. 29 February 1988, 0856 UTC 

b. 29 February 1988, 1707 UTC 

c. 1 March 1988.0833 UTC swath 

Areas of rejections 1 section 4.2 1 are shaded grey. The dashed .me represents 
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Fig. 7 : 1000-850 hPa 31 geopotential thicknesses (in decameter), 29 February 1988, 0856 
UTC. Areas of rejections are shaded grey. 

Fig. 8 : Same as Fig. 7, but 500-300 hPa. 

Fig. 9 : 31 total precipitable water (in kg/m 2 ) 

a. 29 February 1988, 0856 UTC 

b. 1 March 1988, 0833 UTC 
Areas of rejections are shaded grey. 

Fig. 10 : same as 9 but for the layer 1000-850 hPa ; values are muliplied by 10 
Fig. 1 1 : same as 10 but for the layer 850-500 hPa 

Fig. 12 : Composite SSM/I imagery of the cold air outbreak. Integrated water vapor 
derived from Eqn. 1. Values range from 1 to 15 kg/m 2 . The dashed line represents 
approximately the sea ice edge. 

a. 29 February 1988, easternmost swath is from 1624 UTC, westernmost is from 
2124 UTC. 

b. 1 March 1988, easternmost swath is from 0041 UTC, westernmost is from 0542 

UTC. 

c. 1 March 1988, easternmost swath is from 1608 UTC, westernmost is from 2111 

UTC. 

Fig. 13 : same as Fig. 12, but for surface wind speeds, derived from Eq..2 and only for 
evening overpasses of 29 February (a) and 1 March (b). Values range from 0 to 20 m/s. 
The absence of determination ( grey areas over sea ) indicates the presence of precipitation. 
The dashed line represents the approximate sea ice edge. 

Fig. 14 : Cloud liquid water (a) and scattering index by large ice panicles (b) for SSM/I 
pass 3602, 1 March 1988, 0401 UTC. A polar low is identifiable between 72 and 75 N at 
about 17 E. Note that the intense convective clouds close to the Norwegian Sea are well 
delineated by the scattering index. 
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j-ble I Characteristics of TOYS sounding channels. 
































Table 2 : List ofSSM/I swaths (date, orbit number and time UTC ) and corresponding 
significant meteorological features 


date 

orbit 

time 

characteristics 

28 February 

3581 

1633 


(evening) 

3582 

1815 

early stage of the cold air outbreak 


3583 

1955 



3584 

2135 


29 February 

3587 

0233 


(early 

3588 

0413 

cold air outbreak well established 

morning) 

3589 

0555 


29 February 

3595 

1624 


(evening) 

3596 

1805 

strongest temperature and moisture gradients 


3597 

1947 



3598 

2124 


1 March 

3600 

0041 


(early 

3601 

0221 

weakening phase of the outbreak 

morning) 

3602 

0401 

polar low detection by SSM/I scattering index 


3603 

0542 


1 March 

3609 

1608 


(evening) 

3610 

1750 

polar low signature gone 


3611 

1931 

end of the cold air outbreak 


3612 

2111 
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ABSTRACT 

With the atmospheric water-vapor content information available from the SEASAT and Nimbus-7 Scanning 
Multichannel Microwave Radiometers (SMMR), differences in water-vapor distribution between cyclonic storms 
in different regions of the global ocean can be examined in more detail than previously possible from radiosondes. 
SMMR-derived integrated water vapor is a robust and dependable variable of the same accuracy as integrated 
radiosonde soundings. In this study, maximum and minimum water-vapor content in the vicinity of cold fronts 
of 80 storms that occurred in the North Atlantic, North Pacific and Southern oceans are compared. North 
Atlantic storms were found to have significantly higher maximum and minimum water-vapor content near 
cold fronts on average than North Pacific or Southern ocean storms for both the warm and cold seasons. These 
differences are attributed to warmer sea surface temperatures and air temperatures in the North Atlantic, and 
higher baroclinity and consequently stronger upward motion in North Atlantic storms. Additionally, some of 
the differences may be attributed to the fact that the North Atlantic storms generally occur at lower latitudes 
than the storms in the other regions. Furthermore, the North Pacific storms were found to have significantly 
higher maximum and minimum water- vapor content near cold Fronts on average than the Southern Ocean 
storms for both the warm and cold seasons. These differences are attributable to warmer sea surface temperatures 
in the North Pacific during the warm season, and to less moisture transport by Southern Ocean storms and the 
poleward location of the Southern Ocean storms compared to North Pacific storms during the cold season. Two 
examples of water-vapor content in a South Atlantic storm are given to contrast with the Southern Ocean cases. 

The South Atlantic storm had much higher maximum water-vapor content near the cold front than most 
Southern Ocean storms. 


1. Introduction 

The initiation and development of midlatitude cy- 
clones often occur over oceanic regions. Unfortunately, 
the available in situ information regarding the state of 
the atmosphere in vicinity of these storms is sparse, 
and in large regions in the Southern Hemisphere it is 
almost nonexistent. Since the advent of satellite remote 
sensing, this concern has been alleviated. Satellite-de- 
rived fields, such as temperature soundings, have been 
used as input to global numerical weather prediction 
models and have been found to improve forecasts sig- 
nificantly over data sparse regions, especially in the 
Southern Hemisphere (Bourke et al. 1982; Kalnay et 
al. 1985). Visible and infrared satellite pictures are 
routinely used by forecasters and scientists to track 
storm movement and to observe their development 
and decay (Anderson et al. 1969; Weldon 1977; Smi- 
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gielsky and Elrod 1 985 ) , especially over ocean regions. 
Recently, instruments sensitive to the microwave por- 
tion of the electromagnetic spectrum have been placed 
on polar orbiting satellites. Microwave remote sensing 
techniques offer improvements over previous visible 
and infrared techniques by allowing quantitative mea- 
surements of surface and atmospheric variables even 
in the presence of clouds. Atmospheric variables, such 
as boundary layer winds and integrated water vapor, 
have been successfully obtained from satellite-borne 
microwave instruments and proven useful in the re- 
search of midlatitude cyclones (McMurdie and Kat- 
saros 1985; Katsaros and Lewis 1986; and McMurdie 
et al. 1987; Levy and Brown 1986; Brown and Levy 
1986; Levy 1989). 

Atmospheric water vapor plays an important role in 
midlatitude cyclones. The formation of precipitation 
depends on the supply of water vapor ; it is an important 
consideration in determining the structure of atmo- 
spheric boundary layers around frontal zones (Fleagle 
et al. 1988"), and it can influence the development of 
midlatitude cyclones through the release of latent heat 
(Danard 1964; Reed et al. 1988). Detailed structure 
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and magnitude of the atmospheric water- vapor distri- 
bution in midlatitude cyclones were unavailable before 
microwave remote sensing techniques were developed. 
In this paper, the water-vapor information derived from 
the SEASAT and Nimbus-7 Scanning Multichannel 
Microwave Radiometers (SMMR) is examined to 
study the nature of the atmospheric water- vapor avail- 
ability in midlatitude cyclones. The goal of this paper 
is to document typical water-vapor distributions in 
midlatitude cyclones in specific oceanic regions around 
the globe. Differences in water-vapor distribution in 
storms occurring in various regions are highlighted. 

In the following sections, background concerning 
the satellites and the instruments used and the methods 
of the analysis are discussed. The results are highlighted 
with examples of atmospheric water-vapor distribution 
in midlatitude cyclones occurring in the North Atlantic, 
North Pacific, Southern Ocean, and South Atlantic. A 
summary of all the cases examined in these regions is 
discussed in the last section. 

2. Background 

The SEASAT satellite, launched in June 1978, was 
a unique research satellite devoted entirely to the re- 
mote sensing of the world’s oceans. It carried several 
microwave instruments, and demonstrated that the 
state of the ocean environment can be monitored from 
space (Born et al. 1979). Data collected during this 
mission are still being utilized today, despite the pre- 
mature end of SEASAT only three months after launch. 
Soon after the failure of SEASAT, the Nimbus-7 sat- 
ellite was launched in October 1978. The eight instru- 
ments aboard the Nimbus-7 have provided data con- 
tinuously for over eight years. It has allowed an excep- 
tional opportunity to remotely monitor the state of the 
earth-atmosphere system over a long time period. 

Both the SEASAT and the Nimbus-7 satellites car- 
ried a SMMR instrument. It measures the thermal 
emission of the earth-atmosphere system at five mi- 
crowave frequencies (6.6, 10.7, 18.0, 21.0, and 37.0 
GHz) and in both horizontal and vertical polarization 
on a spatial resolution of 50 km X 50 km for the 18.0- 
and 21.0-GHz channels and 27 km X 27 km for the 
37.0-GHz channel. From this information, quantitative 
measurements of sea surface temperature, surface wind 
speed, sea ice concentration, integrated water vapor, 
integrated liquid water, and rainrate can be obtained. 
Of these parameters, the atmospheric water-vapor 
measurements have proven most reliable and as ac- 
curate as radiosonde derived estimates of integrated 
water vapor for a wide range of atmospheric conditions 
(Katsaros et al. 1981; Alishouse 1983; Chang et al. 
1984). 

SMMR-derived integrated water- vapor information 
has been useful for a variety of research objectives. 
Chang et al. (1984) and Prabhakara et al. ( 1982) ex- 
amined the monthly, seasonal, and interannual vari- 


ability of precipitable water over the global oceans. 
They identified several prominent features indicative 
of the large-scale atmospheric circulation, including dry 
regions associated with subtropical anticyclones and a 
region of very high atmospheric water-vapor over the 
western Pacific associated with the Walker circulation. 
Prabhakara et al. (1985) compared the water- vapor 
distribution determined from SMMR during normal 
years to the anomalous El Nino year of 1982/83 and 
found that the water-vapor distribution gave a very 
clear picture of the changes in the atmospheric circu- 
lation patterns during an El Nino year. Francis ( 1987) 
also demonstrated the evolution of the monthly vari- 
ation in the global distribution of water vapor for a 5- 
yr period using recalibrated brightness temperatures of 
Nimbus-7 SMMR. The integrated water-vapor distri- 
bution from SMMR has also been used to examine 
individual midlatitude cyclones over the oceans. The 
fields can be used to accurately locate the position of 
surface frontal zones and frontal waves (McMurdie 
and Katsaros 1985; Katsaros et al. 1989). When the 
SMMR water-vapor field is used in conjunction with 
scatterometer derived boundary-layer winds, qualita- 
tive inference of the relative distribution of water vapor 
with height can be made (McMurdie et al. 1987). 

3. Data analysis 

The SEASAT data were obtained from the NASA 
Jet Propulsion Laboratory (JPL) in the form of geo- 
physical parameters including integrated water vapor. 
The water-vapor parameter was calculated by JPL from 
the brightness temperatures of the 18-, 21-, and 37- 
GHz channels using the Wilheit and Chang (1980) 
algorithm. Some minor corrections to the coefficients 
were made after the data were compared to radiosonde 
soundings (SEASAT, 1981). The Nimbus-7 data were 
obtained from Goddard Space Flight Center in the form 
of brightness temperatures. Corrections outlined by 
Kim et al. ( 1984) were applied to the data to account 
for instrument warm up time, cross-swath biases, and 
drift in the 21 -GHz horizontal channel. The corrected 
brightness temperatures were then processed through 
the Chang et al. ( 1984) water- vapor algorithm. 

With the SEASAT SMMR and Nimbus-7 SMMR 
data, over 80 individual cyclones spanning the period 
September 1978 (SEASAT) and 1979-1982 (Nimbus 
7) are examined. Since the Nimbus-7 SMMR operated 
on an alternate day duty cycle, and after 1 982 the data 
quality deteriorated (see Francis 1987), the number of 
cases available for this study is limited. Individual sat- 
ellite passes that occurred during the same season, and 
in the same oceanic region are grouped together into 
one category. Surface and upper-level analyses from 
the National Meteorological Center (NMC) and the 
Australian Bureau of Meteorology (ABM) and any 
available visible or infrared satellite images are used to 
locate surface fronts. Since the NMC and the ABM 
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surface analyses are used primarily in an operational 
setting, the frontal locations are sometimes inaccurate. 
There are occasions when there are minor discrepancies 
between the SMMR water-vapor patterns and the 
NMC or ABM fro nt l ocation s. When there are dis- 
agreements, they will be mentioned in the subsequent 
analysis. Only satellite passes over cold fronts near the 
surface low of mature cyclones are included. Mature 
cyclones are defined here as storms where the surface 
fronts have occluded and the satellite pictures show a 
distinctive large-scale comma shape. Passes over de- 
caying storms and cold fronts far away from the main 
low pressure center are not included. In addition, cases 
with severe rain interference problems were excluded. 
Rain interference was inferred based on the 37-GHz 
polarization differences and on the water-vapor pat- 
terns themselves. If there was rain interference, it would 
be manifested as low polarization differences in the 37- 
GHz channels and as sporatic anomalously low water- 
vapor amounts within a field of high water- vapor con- 
tent. Very few passes exhibited rain interference, al- 
though rain was present in many cases. A total of seven 
passes were excluded from the study. 

The data were grouped according to season, where 
“warm” season includes storms occurring in Septem- 
ber-November in the Northern Hemisphere and Jan- 
uary-March in the Southern Hemisphere, and the 
“cold” season included storms during December- 
March in the Northern Hemisphere and July-Septem- 
ber in the Southern Hemisphere. During the midsum- 
mer period, June- August^ the storm frequency in the 
North Pacific and North Atlantic is low. We attempted 
to include some storms in this period, but found few 
cases available where the SMMR sampled a storm. Be- 
cause of this, the warm season is slightly different in 
the Southern Hemisphere and the Northern Hemi- 
sphere with respect to the solstice. However, it is likely 
that these results would not be altered substantially if 
we were able to include cases symmetric about the sol- 
stices. This issue is discussed further in section 5. 

Three oceanic regions are studied in detail: the North 
Atlantic, the North Pacific, and the Southern Ocean. 
A few passes from the South Atlantic were also studied, 
but this region is near where the Nimbus-7 SMMR was 
turned on and off and the usable samples were fewer. 
The exact sampling regions are shown in the subse- 
quent sections. In Table 1, the time periods used in 
this study and the number of passes included during 
each period for each area and season are given. A total 
of 1 80 satellite passes are used in the analysis, an av- 
erage of more than two passes per storm. Frontal anal- 
yses shown in all the figures are taken from either NMC 
or ABM unless otherwise noted. 

For each case, a representative maximum and min- 
imum integrated water-vapor content is picked and 
cases in the same category are averaged together. The 
method for selecting the maximum and minimum wa- 
ter-vapor values is as follows. The maximum value is 


either the highest contour in the vicinity of a front (i.e., 
within 200 km of the front) or 2 kg m -2 higher than 
the highest contour if there is a substantial region with 
values higher than the highest contour. The latitude at 
which the maximum value was obtained was as near 
as possible, depending on the sampling by the satellite, 
to the surface low or the triple point, where the oc- 
cluded, warm and cold fronts coincide. The minimum 
value is the value of the lowest contour within 500 km 
of the cold front, towards the cold air, or 2 kg m~ 2 
lower than the lowest contour depending on which is 
more appropriate. 

In addition, statistical tests were made to see if the 
average or mean of one category was significantly dif- 
ferent than the mean of another category using the two 
sample /-test. A minimum confidence level of 95% was 
chosen as the requirement for a difference between two 
categories to be considered significant. However, in 
most of the subsequent analysis, the difference between 
two categories was often significant at the 99% confi- 
dence level. 

Every effort was made to insure that the sampling 
was as random as possible. Every time the satellites 
went over any cold frontal zone and there was no severe 
rain interference or the storm was not decaying, the 
case was included in the statistics. Again, the sampling 
limitations of the Nimbus-7 satellite limited the num- 
ber of cases available for each category. Since several 
years worth of data were included, the sample size was 
sufficient for meeting the significance criteria. 

4. Results 

a . North Atlantic 

Storms that developed in the region east of the east 
coast of the United States from 25° to 55 °N and 80° 
to 20 °W were considered to be North Atlantic storms. 
A map of this region is given in Fig. I. In addition, the 
locations of the warm and cold season storms are in- 
dicated on this map. The western Atlantic Ocean along 
the east coast of the United JStates is a preferred region 
for the development of rapidly deepening cyclones (i.e., 
storms that during their deepening stage had pressure 
falls greater than I mb h~ ! for 24 hours; Sanders and 
Gyakum 1980). Since it is of interest to explore the 
possible differences in SMMR-derived water-vapor in 
these storms compared to regularly developing storms, 
the storms are separated into the categories “rapidly 
deepening” and “standard”. 

A typical example of the water-vapor distribution 
for a warm season North Atlantic mature cyclone is 
given in Fig. 2. The region of maximum water vapor 
is in the vicinity of the cold front and the amounts are 
quite high. There is a strong gradient behind the cold 
front to the northwest of the area of maximum water 
vapor, and the minimum water- vapor values are 15- 
16 kg m' 2 . There is an Increase of water vapor to 22 
kg m -2 in the vicinity of the low center. South of the 
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Table L Dates and number of satellite passes over cold fronts occurring in different regions used in the study.* 


North Atlantic 


North Pacific 


Southern Ocean 


Dates 

Number 

Dates 

Number 

Dates 

Number 

September 1978 

6 

September 1978 

8 

January-March 1979 

27 

November 1980 

12 

November 1980 

10 



September 1981 

4 

September 1979 

5 



November 1981 

l 

September 1981 

3 



Totals 

23 


26 


27 

January-March 1979 

20 

January-March 1979 

18 

September 1978 

9 

January-March 1981 

19 

January 1981 

6 

July 1979 

18 

January 1982 

3 

January, March 1982 

4 

September 1979 

3 

Totals 

42 


28 


30 

Grand total 

176 






* Passes that sampled the warm sector, cold fronts, and the cold sector are all included in this table. This differs from Figs. 1, 4, and 7 
which include the locations of passes sampling cold fronts only. 


cold front, the water-vapor pattern exhibits bands of 
dry and moist air. The band of moist air at 35 °N, 40° W 
is associated with an old, stationary front. The expla- 
nation for the band of high moisture content between 
the cold front and the stationary front is not known at 
this time. 

A typical example of the moisture distribution in a 
cold season mature cyclone in the North Atlantic is 
given in Fig. 3. There is a maximum in water- vapor 
content along the cold front of 32 kg m“ 2 and high 
values of 26 kg m -2 along the warm front. There is 
also a local maximum along the occluded front of 20 
kg m~ 2 . The region of minimum water vapor in the 
cold sector was not sampled. 

A summary of the average maximum and minimum 
water-vapor content of rapidly deepening and standard 
storms are presented in Table 2. The two-sample r-test 
was used to test for significant differences in the max- 
imum and minimum water-vapor contents between 
rapid and standard storms of the same season. These 
tests demonstrated that there is no significant differ- 
ences in maximum water- vapor content in the vicinity 
of mature cold fronts in rapidly deepening cyclones 
and standard cyclones. The difference between the av- 
erage minimum water-vapor content between rapid 
and standard storms during the warm season is signif- 
icant at the 95% confidence level, but the difference is 
not significant during the cold season. However, the 
average maximum and minimum water- vapor contents 
in the warm season storms (total storms category) are 
significantly higher than in cold season storms. This is 
largely a consequence of the difference in the mean air 
temperatures between the seasons. 

Additional passes over cold fronts of developing sys- 
tems and over warm and occluded fronts of developing 
and mature systems were examined by McMurdie 
(1989). In that study, the average maximum water- 
vapor content near cold fronts was significantly larger 
than the average maximum water-vapor content near 


warm fronts, which in turn was larger than the average 
maximum near occluded fronts. This result can be ex- 
plained by the fact that the occluded fronts lie farther 
north than warm fronts and warm fronts lie farther 
north than cold fronts during the mature stage of a 
cyclone. 

b. North Pacific 

The North Pacific category included storms in the 
region from 25° to 60°N and from 120° to I80°W. A 
map of this region is given in Fig. 4 along with the 
locations of the different storms sampled. Examples of 
storms that developed west of 1 80° W are discussed by 
McMurdie (1989) but are not included here. 

A typical example of water- vapor distribution in a 
warm season North Pacific mature cyclone is given in 
Fig. 5. In this case, on 1000 UTC 14 September 1981, 
the maximum water vapor lies along the cold front 
and is greater than 40 kg m“ 2 . Near the warm front in 
the warm sector, the values are also high, greater than 
32 kg m -2 . At 50°N, there is a distinct curvature to 
the contours near the front that parallels a wave along 
the front analyzed by NMC. The driest air sampled in 
this pass lies well west of the cold front at 45°N, where 
the integrated water- vapor content is 16 kg m~ 2 . In 
Fig. 5, three x’s indicate where there is suspected rain 
interference. Despite these small regions of heavy rain, 
the overall water-vapor pattern is retained. 

In contrast. Fig. 6 is a typical example of water-vapor 
distribution during the cold season in the North Pacific. 
Here, again, the maximum water-vapor content lies 
along the cold front; however, the amount is only 28 
kg m’ 2 . The region of minimum water vapor is in a 
broad region immediately behind and northwest of the 
cold front and the amounts are 6-8 kg m -2 . Both of 
these regions have significantly lower water- vapor con- 
tent than similar localities in the warm season case, as 
expected. During the warm season the air temperatures 
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(a) Warm Season 



(b) Cold Season 



Fig. 1. Map of the North Atlantic region. Storms existing between 80° to 20°W and 25° to “ 
60 e N were included in this study, (a) Locations where the satellite sampled the cold fronts of 
rapidly deepening cyclones ( R ) and standard cyclones ( S ) during the warm season, (b) Locations 
where the satellite sampled the cold fronts of rapidly deepening cyclones (R) and standard cyclones 
(S) during the cold season. 


are significantly wanner enabling the atmosphere to 
hold more water vapor, which is clear in these figures. 

The average maximum integrated water- vapor con- 
tent of 28 passes over mature cold fronts during the 
warm season in the North Pacific was 34.9 kg m" 2 and 
the standard deviation was 6.G kg m‘ 2 . The average 


minimum was 1 3.0 kg m~ 2 and the standard deviation 
was 3.2 kg m -2 . The average maximum of 28 passes 
over cold fronts during the cold season was 25.6 kg 
m 2 and the standard deviation was 5.6 kg m“ 2 . The 
average minimum was 9.8 kg m -2 and the standard 
deviation was 3.5 kg m~ 2 . Using the two-sample Mest, 
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are obtained from NMC operational surface analysis. 


the average maximum and minimum water-vapor 
content in the vicinity of cold fronts of mature systems 
is significantly higher during the warm season than 
during the cold season. As with the North Atlantic 
cases, this is largely due to the colder air and sea surface 
temperatures during the cold season. 

Additional passes over developing systems and over 


warm and occluded fronts were examined by Mc- 
Murdie ( 1989 ). In that study, cold fronts were found 
to have significantly higher maximum water-vapor 
content than warm fronts, and warm fronts were found 
to have significantly higher maximum water-vapor 
content than occluded fronts in both seasons. However, 
no significant difference was found between average 
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Fig. 3. An example of SMMR integrated water- vapor distribution 
for 1200 UTC 7 February 1979 in the cold season for the North 
Atlantic. Contours are every 4 kg m' 2 . Frontal positions are obtained 
from NMC operational surface analysis. 


maximum water- vapor content near cold fronts of de- 
veloping cyclones and cold fronts of mature cyclones. 
Interestingly, average maximum water- vapor content 
near warm fronts of developing systems was signifi- 
cantly higher than average maximum water-vapor 
content in warm fronts of mature systems. In this re- 


gion, it is likely that the typical warm front of mature 
systems lies significantly poleward of the typical warm 
front of developing systems, whereas the average lati- 
tude of the sampled cold fronts was similar for devel- 
oping and mature systems. 

c. Southern Ocean 

Storms that occurred in the region 90°-160°E and 
30°-60°S were considered to be Southern Ocean 
storms. This area includes the ocean west and south 
of Australia, as seen in the map given in Fig. 7. The 
locations of the sampled storms in this study are also 
indicated on the map. The operational analyses of the 
ABM were used to determine the location of the fronts. 

A typical example of a warm season cyclone in the 
Southern Ocean is given in Fig. 8. This Nimbus-7 ov- 
erpass occurred at 0521 UTC 14 January 1979 just 
west of the west coast of Australia. The region of max- 
imum water vapor lies in a relatively narrow band along 
the cold front with values of 24 kg m“ 2 . At the cold 
front near 50°S the maximum water-vapor content is 
less at 16-20 kg m -2 . Immediately west of the cold 
front the water vapor decreases rapidly to a minimum 
of 11-12 kg m" 2 . Near the low pressure center, the 
water- vapor content increases to 16-17 kg m -2 . Pole- 
ward of the system at 55 °S, the water- vapor content is 
very low, less than 6 kg nT 2 . The region south of 55°S 
is often very cold and dry in all seasons. 

A typical example of a cold season Southern Ocean 
cyclone is given in Fig. 9. Again, the maximum water 
vapor lies in the vicinity of the cold front, however, 
the amounts are very low, 20 kg m -2 . There is a fairly 
strong decrease behind the cold front, to the south, to 
a region of minimum water-vapor content of 4-5 kg 
rrT 2 . Although the water- vapor pattern is typical of 
midlatitude storms in many regions, the water-vapor 
content is exceptionally low. 

The average water-vapor content for warm and cold 
season storms are summarized as follows: of 27 satellite 
passes over warm season cyclones, the average maxi- 
mum in the vicinity of cold fronts is 25.3 kg in' 2 and 


Table 2. Average maximum and minimum water-vapor content for warm and cold season cold fronts of North Atlantic storms. 
Average water-vapor content and standard deviation are given in kilograms per square meter (kg m -2 ), and n is the number of cases. 




Maximum 



Minimum 


Case 

Average 

Standard 

deviation 

n 

Average 

Standard 

deviation 

n 

Warm season 
Standard 

41.5 

4.5 

12 

18.1 

4.5 

13 

Rapidly deepening 

43.0 

4.7 

6 

14.8 

4.4 

10 

Total 

42.0 

4.5 

18 

16.6 

' 4.6 

23 

Cold seaspjf: ^ - 
Standard - 

30.6 

5.2 

27 

14.5 

4.7 

29 

Rapidly deepening 

33.5 

4.7 

11 

14.1 

5.6 

10 

Total 

3L5 

5.0 

38 

14.3 

4.9 

39 
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(a) Warm Season 



(b) Cold Season 



Fig. 4. Map of the North Pacific region. Storms existing between 180° to 120°W and 25° to- 
60°N were included in this study. Locations where the satellite sampled the cold fronts are indicated 
by a star (*) during (a) the warm season and (b) the cold season. 


the standard deviation is 3.9 kg m“ 2 . The average min- 
imum behind cold fronts is 9.0 kg m" 2 and the standard 
deviation is 1 .9 kg m" 2 . Of 30 satellite passes over cold 
season cyclones, the average maximum in vicinity of 
cold front is 20. 1 kg m“ 2 , the standard deviation is 4.5 
kg m -2 . The average minimum behind cold fronts is 
8,2 and the standard deviation is 2.7. As in the other 
regions, there is a significant difference in the maximum 


water-vapor contents between the cold season and 
warm season storms. However, the difference is not 
nearly as large as in the North Pacific or in the North 
Atlantic cases. 

d. South Atlantic 

The Southern Ocean cases presented in the previous 
section are not representative of the entire Southern 
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Fig. 5. An example of SMMR integrated water- vapor distribution 
for 1000 UTC 14 September 1981 in the warm season for the North 
Pacific. Contours are every 4 kg m' 2 . Frontal position is obtained 
from NMC operational surface analysis. The small x’s indicate regions 
of suspected rain interference. 


Hemisphere oceanic regions. In this section two ex- 
amples from a South Atlantic ocean storm are pre- 
sented to contrast with the Southern Ocean cases. Rig- 
orous tests regarding significant differences between 
storms in the Southern Ocean and this region have not 
been performed. 

In Figs. 10 and 1 1, a storm that developed off the 
east coast of South America is sampled twice by Nim- 
bus-7 SMMR, two days apart. In the earlier pass (Fig. 
10), 1400 UTC 28 January 1979, there is a large 


amount of water vapor in the vicinity of the cold front, 
approximately 32 kg m~ 2 at 45 °S, and the water vapor 
increases equatorward to greater than 40 kg m“ 2 . There 
is no distinct maximum along the front. The gradient 
behind the front, towards the south, is strong, and the 
minimum water-vapor content in this region is 1 1 kg 
m“ 2 . In the later pass (Fig. 1 1 ), 1 300 UTC 30 January 
1979, the vapor content in the vicinity of the cold front 
is higher, greater than 40 kg m“ 2 at 40 °S, and the 
amount decreases slightly equatorward. There is a very 
strong gradient poleward to a minimum amount of 7 
kg m -2 . The front position was taken from the 1200 
UTC 30 January 1979 ABM map. It appears that the 
ABM frontal location is too far into the cold dry air. 
Our past experience would place the front closer to the 
region of maximum water- vapor content along the 24 
or 28 kg m~ 2 contour line in this figure. 

In both of these cases the maximum water-vapor 
content is much higher than in any of the Southern 
Ocean storms. On seasonal time scales, the region east 



FlG. 6. An example of SMMR integrated water- vapor distribution 
for 2200 UTC 15 February 1979 in the cold season for the North 
Pacific. Contours are every 4 kg m~ 2 . Frontal positions are obtained 
from the NMC operational surface analysis. 
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(a) Warm Season 



(b) Cold Season 


FIG. 7. Map of the Southern Ocean region. Storms existing between 90° to 160°E and 25° to 
60 S were included in this study. Locations where the satellite sampled the cold fronts are indicated 
by a star (* ) during (a) the warm season and (b) the cold season. 



of South America and equatorward of 40°S has total 
water-vapor content greater than 30 kg m" 2 during De- 
cember-February ( see Fig. 12). This is higher than the 
vapor content at the same latitude in the vicinity of 
Australia (15-20 kg m~ 2 ). In addition, this particular 
area is a known region of cyclogenesis (Physick 1981; 
Carleton 1980). Using seasonal mean fields for the pe- 
riod July-August 1982, James and Anderson (1984) 
found a high incidence of transient eddy activity 


(baroclinic waves) over the South Atlantic, east of 
South America. They found that this activity could not 
be explained by low- and midlevel baroclinicity defined 
on the basis of dry, linear theory. They suggested that 
an increase of moisture entrainment into the westerlies 
in this region could account for the large transient ac- 
tivity there. They explored this possibility using model 
derived moisture, but lacked sufficient moisture ob- 
servations to substantiate their results. The SMMR data 
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FiG. 8. An example of SMMR integrated water- vapor distribution 
at 0600 UTC 1 4 January 1 979 in the warm season over the Southern 
Ocean. Contours are every 4 kg m' 2 . Frontal position is estimated 
from the 0000 UTC and 1200 UTC 14 January 1979 ABM opera- 
tional surface analyses. 


support their suggestion and show there is higher 
amount of moisture in storms observed east of South 
America compared to Southern Ocean storms. 

5. Discussion 

In the previous sections, the typical distribution of 
water vapor in midlatitude cyclones in the North Pa- 
cific, North Atlantic, and Southern Ocean were pre- 
sented and two passes over a storm in the South At- 
lantic were discussed. In Table 3, the average maximum 
and minimum water-vapor content in mature cold 
fronts for these three areas are repeated. Using the^®K; 
sample /-test, the North Atlantic cold fronts have sig- 
nificantly higher average maximum and minimum 
water-vapor contents than the North Pacific ones for 
both seasons; and the North Pacific cold fronts in turn 


have significantly higher average maximum and min- 
imum water-vapor content than the Southern Ocean 
cold fronts for both warm and cold seasons. 

It is clear from Table 3 that the North Atlantic region 
has very high moisture content in midlatitude storms. 
This can partially be explained by the relatively high 
sea surface temperatures in the North Atlantic and 
consequently higher evaporation and moisture flux into 
the atmosphere in this region, especially in the warm 
season. Consequently, the seasonal mean integrated 
water-vapor content in the North Atlantic is higher 
than in the North Pacific and the Southern Ocean for 
both seasons. During the warm season, it ranges from 
35-20 kg m~ 2 in the latitude interval 30°-50°N in the 
North Atlantic. In contrast, in the North Pacific and 
Southern Ocean, it ranges from 30- 1 5 kg m” 2 and 2 1 - 
1 5 kg m“ 2 , respectively, in the same latitude band (see 
Fig. 12). During the cold season, the mean water- vapor 



Fig. 9. An example of SMMR integrated water- vapor distribution 
at 1600 UTC 3 September 1979 in the cold season over the Southern 
Ocean. Contours are every 4 kg m" 2 . Frontal position is estimated 
from the 1200 UTC 3 September 1979 and the 0000 UTC 4 Sep- 
tember 1979 ABM operational surface analyses. 
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28 JAN 1979 14 GMT 


Fig. 10. Integrated water- vapor distribution for a South Atlantic 
storm that occurred at 1400 UTC 28 January 1979. Contours are 
every 4 kg m" 2 . Frontal position is from the 1200 UTC 28 January 
1979 ABM frontal analysis. 


content in the 30°-50° latitude band ranges from 25- 
15 kg m“ 2 in the North Atlantic, 22-12 kg m“ 2 in the 
North Pacific, and 17-10 kg m“ 2 in the Southern Ocean 
(see Fig. 12). In addition, there is a western boundary 
current in the North Atlantic in the Gulf Stream, which 
is a permanent feature just off the east coast of the 
United States where warm ocean waters travel rapidly 
northward and eastward from the tropics to the mid- 
latitudes. The availability of atmospheric moisture 
from this warm ocean source must contribute to the 
moisture content in weather systems in the North At- 
lantic region. Furthermore, the North Atlantic, specif- 
ically along the East Coast, is a known region of cy- 
clogenesis; often the storms that form there deepen 
rapidly. The proximity of the warm Gulf Stream to a 
source of cold air of continental origin enhances low- 
level baroclinicity and decrease the static stability in 
the atmosphere producing a favorable environment for 


cyclogenesis (Hadlock and Kreitzberg 1988). Since it 
is such an active region of cyclogenesis and frontogen- 
esis, the synoptic and mesoscale convergence and up- 
ward vertical motion in the vicinity of fronts may be 
stronger than in most storms, thus, creating higher in- 
tegrated moisture in the frontal zones of these storms. 
In addition, the North Atlantic storms used in this study 
occurred at lower latitudes than the North Pacific and 
Southern Ocean storms. When the latitude of the max- 
imum water-vapor content in the vicinity of the cold 
front for each pass over a storm is used to indicate the 
location of each storm, the average latitude of the North 
Atlantic storms is 33°N. In contrast, the average lati- 
tude of the North Pacific and Southern Ocean storms 
is 38°N and 42 °S, respectively. All these factors must 
contribute to the high moisture content in storms in 
the North Atlantic. 

The integrated moisture content in the region of the 
minimum behind the cold front was also highest in 
the North Atlantic for both seasons. Although the air 
in the region where the water vapor is minimum typ- 
ically originates from very cold and dry continental 



30 JAN 1979 13 GMT 


Fig. 1 1. Integrated water-vapor distribution for a South Atlantic 
storm that occurred at 1300 UTC 30 January 1979. Contours are 
every 4 kg m" 2 . Frontal position is from the 1200 UTC 30 January 
1979 ABM frontal analysis. 
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areas, when it pours over the relatively warm ocean, 
the surface flux of heat and moisture modifies this air 
strongly. Petterssen et al. ( 1 962 ) calculated the sum of 
sensible and latent heat flux in the region west of a 
typical North Atlantic storm and found that it ranged 
from 100-400 W m“ 2 . In addition, Raman and Rior- 
dan ( 1988 ) found during the Genesis of Atlantic Lows 
Experiment (GALE) that a sum of sensible and latent 
heat flux of 800 W m" 2 was common during cold-air 
outbreaks following midlatitude storm passages. This 
vigorous exchange of heat and moisture from the ocean 


to jhe cold air behind fronts is probably reflected in 
the slightly greater integrated water-vapor content in 
the minimum region in North Atlantic storms. 

As mentioned in section 3, the sampling periods 
representing the Southern Hemisphere warm season 
and the Northern Hemisphere warm season do not 
have identical separation from the solstices. This was 
due to a lack of cases av ail able in the Northern Hemi- 
sphere during the early warm season (July and August): 
However, if there were enough storms sampled during 
this time period, one would expect the maximum wa- 
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Table 3. Average maximum and minimum water-vapor content for cold fronts in mature cyclones. 

Average water-vapor content and standard deviation are given in kilograms per square meter (kg m 2 ), and n is the number of cases. 


Maximum Minimum 


Case 


Standard 

Average deviation 


Standard 

n Average deviation n 


Warm season 


North Atlantic 

42.0 

4.5 

North Pacific 

34.9 

6.0 

Southern Ocean 

25.3 

3.9 

Did season 



North Atlantic 

31.1 

5.0 

North Pacific 

25.6 

5.6 

Southern Ocean 

20.1 

4.5 


18 

16.6 

4.6 

23 

28 

13.0 

3.2 

24 

27 

9.0 

1.9 

26 

38 

14.3 

4.9 

39 

27 

9.8 

3.5 

28 

30 

8.2 

2.7 

29 


ter-vapor content in the vicinity of the cold fronts of 
early warm season storms to be the same as or greater 
than the water- vapor content in cold fronts of late warm 
season storms examined here. This is partly because 
the average water-vapor content in the North Pacific 
and North Atlantic is higher during the early warm 
season months (June-August) than during the late 
warm season months (September-November) (e.g.. 
Fig. 2 of Prabhakara et al. 1985). On the other hand, 
in the Southern Ocean, the average water- vapor content 
does not change appreciably from early to late warm 
season and the storm frequency and intensity is about 
the same during the early and late warm season. 
Therefore, one would expect that the average waters 
vapor content in the vicinity of cold fronts in the 
Southern Ocean calculated in this study to be repre- 
sentative for the entire warm season. If the fortunate 
circumstances had existed so that the Southern Hemi- 
sphere and Northern Hemisphere warm seasons sam- 
ples had been so chosen as to be exactly equidistant 
from the respective solstices, we have confidence that 
our general conclusions would have been similar. 

In both the warm and cold seasons, the Southern 
Ocean average maximum and minimum water- vapor 
contents stand out as exceptionally low compared to 
the North Pacific and North Atlantic cases. In some 
respects, this could be due to a different large-scale cir- 
culation in the Southern Hemisphere. Several authors 
using a variety of data sources and analysis methods 
have arrived at the following characteristics of the gen- 
eral circulation (Trenberth 1980; Physick 198.1 ; James 
and Anderson 1984): l) the seasonal patterns of the 
surface pressure and 500-mb height are dominated by 
wavenumber 3 configuration with troughs typically at 
longitudes 1 10°E, 120°W, and 40° W; 2) the seasonal 
cycle is much less pronounced than in the Northern 
Hemisphere with a high frequency of transient cyclones 
throughout the year; and 3) the region of low-level 
baroclinicity and the high-latitude wintertime jet is 
further poleward than in the Northern Hemisphere. 
The most important feature that contributes to inte- 


grated water-vapor content of these three is the last, 
which implies that the jet and the storm and depression 
tracks tend to be at high latitudes in the Southern 
Ocean. This would presumably reduce the transport 
of warm moist air from lower latitudes by storm sys- 
tems. 

During the warm season, the Southern Ocean does 
not warm substantially. At 40 °S it is roughly 1 5°C and 
at 50°S it is 10°C; whereas the North Pacific is 20°C 
and 14°C at these same latitude bands (van Loon 
1984). In addition, the warm season integrated water- 
vapor content in the Southern Ocean ranges from 2 1 
kg m" 2 at 30°S to 15 kg m" 2 at 50°S; whereas in the 
North Pacific, the range is from 30 kg m~ 2 at 30°N to 
15 kg m“ 2 at 50°N (see Fig. 12). It appears that the 
lower sea surface temperatures and therefore lower at- 
mospheric temperatures and mean atmospheric water- 
vapor contents in the Southern Ocean summer is 
probably the most important reason why there is such 
low moisture content in the Southern Ocean storms 
during the warm season. 

It is somewhat puzzling that the Southern Ocean 
winter storms have significantly less total moisture than 
the North Pacific storms. In both regions the storms 
are predominately in the mature stages of development 
and the sea surface temperatures are climatologically 
the same, ranging from 17°C at 35° latitude to 7°C at 
50° latitude (van Loon 1984) in both oceans during 
winter. There are several possibilities to account for 
this difference, including: 1) storms in the Southern 
Ocean occur significantly poleward of their North Pa- 
cific counterpart; and 2) Southern Ocean storms trans- 
port significantly less moisture poleward than the North 
Pacific storms. 

In order to test the hypothesis that the Southern 
Ocean storms occurred significantly poleward of the 
North Pacific storms, the average latitude of the surface 
low pressure center for the Southern Ocean and the 
North Pacific storms used in this study was calculated. 
The average latitude of the Southern Ocean low pres- 
sure centers was 51°S and the average latitude of the 
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North Pacific low pressure centers was 49°N. However, 
the standard deviations were so large ( a = 7.0 for the 
Southern Ocean and a - 6.5 for the North Pacific) 
that the difference between the averages was not sig- 
nificant. In addition, the average latitude of the cold 
fronts sampledln the Southern Ocean was 42°$ and 
the average latitude of the cold fronts sampled in the 
North Pacific cold fronts was 38 °N. Despite the large 
standard deviations (a = 5.5 for the Southern Ocean 
and a = 7.3 for the North Pacific), the averages are 
significantly different at the 97% confidence level. 
Therefore, there is limited evidence that the Southern 
Ocean cold fronts in this study were slightly poleward 
of the North Pacific cold fronts. However, the evidence 
does not appear sufficient to fully account for the dif- 
ferences between the water- vapor content in Southern 
Ocean storms and in North Pacific storms. 

The second hypothesis that the Southern Ocean cold 
season storms transport significantly less moisture than 
the North Pacific storms is more difficult to test Ad- 
equate three-dimensional wind and moisture infor- 
mation for all the cases examined here is lacking. 
However, several authors have computed the integrated 
water- vapor distribution and vertically integrated me- 
ridional and zonal components of moisture flux on 
seasonal time scales. With radiosonde and model de- 
rived information, Howarth ( 1 983 ) calculated that the 
meridional component of moisture flux in the Southern 
Ocean region during winter months was 25 kg m' 1 
s _i . Peixoto et al. ( 1981 ) performed a similar calcu- 
lation using radiosonde derived information and com- 
puted a value of 89 kg m _i s" 1 for the North Pacific. 
They both attributed the poleward meridional flux of 
moisture to cyclone activity. In addition, the semiper- 
manent Aleutian Low may play a role in the higher 
vapor content of North Pacific storms. Rasmusson 
(1972) showed that the zonally averaged meridional 
transport of moisture due to stationary eddies is sig- 
nificant. During the winter season in the Northern 
Hemisphere it accounts for approximately one-third 
of the total meridional transport of moisture. The 
Southern Ocean, on the other hand, does not have 
semipermanent features such as the Aleutian Low. This 
moisture transport may explain why the seasonal total 
integrated water-vapor content is slightly higher in the 
North Pacific basin than in the Southern Ocein fsee 
Fig. 12), especially in the eastern North Pacific. The 
North Pacific storms, therefore, may occur in an en- 
vironment with higher vapor content than the Southern 
Ocean storms on average. In any case, using the values 
calculated by Peixoto et al. (1981) and Rasmusson 
( 1972), the relative contribution to the total meridional 
moisture flux by transient eddies in the North Pacific 
is approximately 60 kg m _1 s“\ twice that calculated 
by Howarth for the Southern Ocean. These studies and 
others, support the hypothesis that Southern Ocean 
cold season cyclones transport less water vapor than 


the North Pacific cold season cyclones. Yet, without 
information on temporal and spatial scales of the cy- 
clones themselves, this must remain an open question. 
Perhaps in the future, with high resolution wind in- 
formation, much like that of the SEASAT Scatter- 
ometer, and high resolution water- vapor information, 
like that of the SMMR or the current Special Sensor 
Microwave Imager (SSM/I), the meridional flux of 
boundary-layer water vapor poleward in cold season 
midlatitude cyclones can be computed for these re- 
gions. 

McMurdie ( 1989) found some interesting differences 
in maximum water-vapor content in cold fronts and 
warm fronts. The North Pacific and North Atlantic 
mature storms were found to have significantly higher 
maximum water- vapor content in cold fronts than in 
warm fronts. In addition, in the North Pacific cases, 
developing warm fronts have significantly more water- 
vapor than mature warm fronts. These observations 
can be partially explained by differences in the mean 
temperature of the air column at mature and devel- 
oping warm and cold fronts. Mature warm fronts tend 
to lie significantly poleward of the mature cold front 
of the same system and of the position of the warm 
front during the storm’s developing stage. In many 
cases, the air at the mature warm front has been mod- 
ified and has cooled over time and mixed with the 
cooler air to the northeast resulting in lower water- 
vapor content at the mature warm front. Observations 
and modeling of the dynamical development of warm 
and cold fronts over time are consistent with the above 
observations ( Eliassen 1962V Hoskins and Heckley 
1981;KeyserandPecnick 1987; and Levy 1989). These 
authors show that cold fronts continue to develop sharp 
temperature gradients and have strong convergence 
(ancFtherefore higher moisture convergence) as the 
storm system develops, whereas warm fronts do not. 

In this study, significant differences in average max- 
imum and minimum water-vapor content were found 
between storms that occur in the North Atlantic, North 
Pacific, and Southern Ocean. In addition, there is lim- 
ited evidence that the South Atlantic storms have much 
higher maximum water- vapor content than the South- 
ern Ocean storms, lending observational evidence in 
support of the speculation of James and Anderson 
(1984) concerning the high transient activity in the 
South Atlantic. These differences in water- vapor con- 
tent are attributed to differences in ocean temperatures 
and moisture advection and convergence. The relative 
importance of these mechanisms in explaining the dif- 
ferences were discussed. These observational results 
comprise the first documentation of typical water-vapor 
contents in midlatitude cyclones occurring in these re- 
gions. Full explanation of the regional and seasonal 
differences in water-vapor content of cyclones must 
await further studies probably incorporating diagnostic 
analysis of limited area or global numerical models. 
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Despite the limitations inherent in the SEASAT and 
Nimbus-7 SMMR, this study demonstrates what can 
be done with the SSM/I data on current and future 
DMSP satellites and eventually with data from the mi- 
crowave radiometers in the Earth Observing System 
planned for 1995 and beyond. 
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ABSTRACT 

The cloud cover and cloud water content obtained from integration of a limited area mesoscale 
prediction model are compared with the corresponding quantities derived from AVHRR and 
SSM/I satellite data. The model has an elaborate condensation-cloud parameterization scheme, 
in which cloud water is a prognostic variable. From AVHRR data, both cloud cover and liquid 
water content are retrieved by utilization of radiances from the visible, near infrared and infrared 
channels. From the 37GHz polarization difference in the microwave data obtained from 
SSM/I measurements, liquid water content is interpreted; this is only possible over the sea. An 
important feature of this comparison exercise is that the retrieved data are distributed on the 
same gnd as the model. The difference between the simulation and the analysed cloud fields can 
thus easily be quantified and exhibited on the model domain. It is demonstrated that this type 
of retrieved meteorological parameters constitutes potentially useful data samples for quan- 
titative verification of model simulations. There are some weaknesses in the retrieval methods 
that have to be removed before we can have sufficient confidence in the data to make quan- 
titative verification meaningful. Furthermore, calibration against cloud parameters obtained 
through other measurements is necessary, although the three data sets show a remarkably good 
agreement in pattern and range of magnitude in the liquid water content. Even if improvements 
are needed, the present study has made it possible to make inferences about some deficiencies 
in the parameterization scheme of the model. 


1. Introduction 

During the past decade, a great number of 
numerical model studies have been devoted to 
investigation of sensitivity to cloud-radiation con- 
ditions in the atmosphere. By and large, it can be 
said that all those studies indicate that there is an 
important impact from the cloudiness on the 
resulting circulation in NWP models and GCM 
(see list of acronyms in Appendix) as well as in 
climate models (Manabe and StoufFer, 1980; 
Ramanathan et al., 1983; Rockner and Schlese, 
1987; Wetheraid and Manabe, 1988; Cess et al., 


* Contribution number 625. 


1989; Slingo, 1990). It can furthermore be stated 
that very little attention generally is paid to the 
model clouds themselves. With regard to the sen- 
sitivity to clouds, it is urgent that much more work 
be done on the problem of condensation-cloud 
parameterization. The more realistic the model 
cloudiness becomes, the more confidently we can 
draw conclusions from model results. In the most 
recent papers, a concern may be found about the 
realism in the model cloud fields (Wetheraid and 
Manabe, 1988; Cess et al, 1989). 

There are two aspects that must be considered in 
a parameterization problem. One is the derivation 
and development of parameterization schemes, the 
other is the verification that model results exhibit 
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a realistic behavior and distribution of, in this case, 
cloudiness. 

The conventional ground-based cloud observa- 
tions are incomplete and spatially too sparse to 
constitute a useful data set in the context discussed 
here. It is not until recently that relevant observa- 
tional data have become available through satellite 
measurements. From the AVHRR data, it is now 
possible to make quantified interpretations of 
cloud cover and cloud water content (Raustein, 
1989). From the microwave measurements, quan- 
tified cloud water content and precipitating water 
content may be deduced (Staelin et al., 1976; 
McMurdie and Katsaros, 1985; Takeda and Liu, 
1987; Katsaros et al., 1989; Petty, 1990; Curry 
et al, 1990; Katsaros et al., 1990). 

The satellite measurements also have their 
limitations. Cloud classification from AVHRR 
data utilizes channels both in the infrared (IR), 
and in the visible part of the spectrum. The inter- 
pretations are sensitive to sun elevation, and 
criteria for distinguishing between clouds and 
snow are not currently included in the method. 
Furthermore, microwave data can only be used 
over ocean areas. 

The interpretation of both AVHRR data and 
microwave data are in an early stage of develop- 
ment. This means that the “observations" obtained 
this way have some degree of uncertainty, espe- 
cially in absolute values. 

There is another important aspect that is com- 
mon to all of these satellite measurements. 
Namely, they give the vertically integrated value of 
the measured quantities. These do not necessarily 
reflect a unique vertical distribution (Slingo and 
Slingo, 1988). 

Because of the development stage of the inter- 
pretation schemes and because of the two-dimen- 
sionality of the measurements, it is important that 
we utilize as many categories of measurements 
as possible for the model verifications. A careful 
comparison of corresponding measurements by 
different methods, as well as comparison of 
model results and measurements, may assist in 
calibrating the observations. In this way an 
enhanced confidence in the observations may 
gradually be achieved. 

In order to verify simulated vertical distribu- 
tions, it hence appears necessary to obtain obser- 
vations also in other ways, such as from dedicated 
field programs such as, for example, FIRE 
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(Cox et al., 1987). Then, as condensation-cloud 
parameterization schemes begin to show an 
acceptable performance, those schemes can be 
utilized to infer vertical distributions from satellite 
measurements. 

There is a long way to go before we have reached 
a satisfactory state with regard to observational 
data of condensation-cloud quantities, as well as a 
realistic simulation of those quantities in various 
models. Considering the attention the cloud topic 
has gained recently, it is indeed urgent to intensify 
the research efforts in both modelling of condensa- 
tion-cloud processes and in the verification of 
those simulations. 

We are therefore of the opinion that, at the pre- 
sent state of the art, cloud parameters, derived 
from satellite data, are not yet quite ready to be 
used for direct absolute verification of model 
simulations. But to improve this situation, it is 
necessary to go through a development phase, in 
which we make comparisons between interpreted 
variables and model variables. This is the purpose 
of the present paper. A study along similar 
lines has recently been published by Saunders 
(1989). 

The simulated cloudiness is produced in a NWP 
model (the AMUB model, which is a version of the 
operational model of the Norwegian Weather 
Service, modified at the University of Bergen), 
which has an elaborate condensation-cloud 
parameterization scheme. 

The observed cloud parameters are obtained 
from two sources. One comprises cloud cover 
and vertically integrated liquid water content, 
objectively derived from satellite AVHRR data 
according to the approach suggested by Raustein 
(1989). The other data set is vertically integrated 
liquid water content, which is obtained by inter- 
pretation of satellite microwave measurements 
according to Petty (1990). 

One model field which would have been 
desirable to compare to a corresponding satellite 
data field, is that of total outgoing longwave radia- 
tion. However, this field is not available from the 
satellite data we are using, as they are taken only 
from a few window bands. 

Specifically, comparisons are made between 
clouds, predicted with the AMUB model and 
satellite clouds at times differing by not more than 
half an hour to one hour. An important feature of 
this study is that the cloud parameters from the 
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satellite data are projected on the same grid as the 
model. 

2. Model features and the derivation of 
cloudiness from AVHRR and SSM/I 
microwave data 

2.1 . Model features 

The AMUB model is a modified version of the 
operational fine mesh model at the Norwegian 
Weather service (DNMI). A description of the 
original version may be found in Gronas et al. 
(1987) and Nordeng (1986). The modifications 
introduced at AMUB are described by Sundqvist 
et al. (1989). Consequently, here we only present 
the broad features of the AMUB model. This 
limited area model is based on the primitive equa- 
tions with sigma as vertical coordinate. In the 
experiments presented here, the horizontal resolu- 
tion is 50 km and there are 10 levels in the vertical, 
with the top at 200 mb. The number of grid points 
at each level is 49 x 61; the geographical area that 
is covered may be seen in Fig. 4. 

The time integration is explicit, but the scheme 
contains spatial and temporal smoothing of the 
tendencies in the staggered grid so that the time 
step is nearly four times as long as the one basically 
obtained from the CFL condition (Bratseth, 1983). 
Vertical eddy diffusion is applied through the 
whole depth of the model atmosphere. The 
exchange coefficient is a function of the 
Richardson number and cloud cover. The sea sur- 
face temperature is prescribed, while over land, an 
energy exchange equation is applied as well as 
a soil moisture formulation, accounting for 
precipitation, storage and run off. Radiative fluxes 
are calculated at all moderievels. The fluxes 
are modified with regard to the current cloud 
cover obtained from the condensation-cloud 
parameterization. The emissivities are not func- 
tions of the cloud water content, but prescribed 
with different values for high, medium high and 
low clouds. 

As indicated, the particular feature of the 
AMUB model is that an elaborate condensation- 
cloud parameterization scheme is employed. This 
scheme carries cloud water as a prognostic 
variable and treats release of latent heat, 
appearance of fractional cloud cover, cloud water 
content and precipitation/evaporation in a consis- 
tent way (Sundqvist, 1988; Sundqvist et al., 1989). 


The parameterization of the microphysical 
processes accounts for coalescence as well as 
Bergeron-Findeisen effects. This scheme is common 
for convective and stratiform condensation, but the 
associated parameters have different values in the 
two stability regimes. 

For the parameterization of convective con- 
densation , a modified Kuo scheme is employed. The 
scheme includes calculation of cloud water. The 
partitioning of converging vapour between moist- 
ening and condensation is essentially governed 
by the relative humidity. The fractional cloud 
cover is a function of the intensity and the depth of 
the convection. The condition for stratiform con- 
densation to appear is that the relative humidity is 
greater than a prescribed threshold value, which is 
different over sea and land. The fractional cloud 
cover is a linear function of the square root of the 
sub-saturation. The cloud wafer advection, which is 
only horizontal, is calculated with an upstream 
scheme to assure positive definite cloud water 
values. 


2.2. Derivation of cloud parameters from A VHRR 
data 

The method of deriving fractional cloud cover 
and vertically integrated liquid water content 
(LWC) is described hy Raustein (1989). A brief 
description of the method is given below. 

Three of the AVHRR channels are utilized, 
viz. channel 1 (0.58-0.68 /im, visible), channel 2 
(0.7-1. 1 pm, near infrared), and channel 4 
(10.5-11.5 /im, infrared). First a statistical method 
called clustering is applied to the data of channels 
1 and 4 to obtain an objective cloud classification. 
In the two-dimensional radiance space, we find 
concentrations of pixels (called clusters), which we 
identify as certain cloud types, or as either land or 
sea surface. Pixels belonging to a certain cluster 
contain clouds of the same type, and the area 
covered by such a pixel is assumed to be overcast 
Pixels that are scattered and thus not belonging to 
a cluster are overcast or partially cloudy. The 
cloud fraction of such a pixel is found in two steps: 
(a) The cloud type, to which the clouds in the pixel 
belong. Is obtained by finding the cloud cluster, to 
which the pixel is closest in the 3-dimensional 
radiance space, where the 3 dimensions are the 
radiances of channels 1, 2 and 4. (b)The cloud 
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fraction is then found by linearly interpolating the 
pixel’s channel 2 radiance between the surface 
cluster’s mean, and the closest cloud cluster’s. 

To compute the LWC and cloud top tem- 
perature we use radiative transfer models of 
Coakley and Chylek (1975) and Arking and Childs 
(1985), together with parameterized LWC and 
channel 4 emissivity according to Stephens (1978). 

In order to have the cloud parameters derived 
from the satellite data projected on the same grid 
as the model, a navigation method developed by 
Kloster at NRSC (personal communication) is 
used to position each pixel in the model grid. The 
grid box value of each cloud parameter is then 
taken as the average of the values of the pixels 
which belong to the grid box. 

2.3. Derivation of liquid water content from 
microwave data 

The microwave data are from the SSM/I on the 
F8 satellite in the US Defense Meteorological 
Satellite Program (DMSP), which was launched 
into a polar orbit in June 1987. This microwave 
radiometer scans a 1390 km wide swath in a 
conical sweep at 53° incidence angle. Operating 
frequencies are 19.35, 22.235, 37 and 85.5 GHz, in 
both horizontal and vertical polarization, except 
the 22.235 GHz channel, which only operates in 
horizontal polarization. The pixel resolution is 
nominally 50, 50, 25 and 15 km, respectively. 
Several possible algorithms for determining cloud 
liquid water content have been developed (e.g., 
Petty, 1990). For the purpose of this model- 
satellite data comparison, a 50-km resolution- 
integrated cloud liquid water algorithm based on 
the 37 GHz polarization was selected. For higher 
spatial resolution and better discrimination, an 
algorithm employing the 85.5 GHz channels has 
been developed (Petty, 1990). Unfortunately, this 
algorithm cannot be used for our August 1988 case, 
due to gain problems with the 85.5 GHz channel in 
vertical polarization, which had started shortly 
before this time. 

The 37 GHz algorithm is based on the theoreti- 
cal depolarization effect of the cloud droplets on 
the polarized emission from the sea surface. A 
relation between the observed polarized signal, P, 
and vertically integrated liquid water mass, L 
( = LWC), is of the following form (Petty, 1990): 

< 2| > 


where tj is effective, microwave transmittance 
associated with the cloud liquid water, a has the 
approximate value 2, k is effective mass extinction 
coefficient of the liquid water in the cloud and 6 is 
the incidence angle (53°). 

The relation for L 37 in terms of polarization at 
37 GHz (Tjtv — Tym) has been found from a 
theoretical brightness temperature model, which 
provides values for a • k: 

L 37 =- 1.42 In(P 37 ) (kgm- 2 ). (2.2) 

P 37 includes corrections for changes in the back- 
ground, that is the wind roughening of the sea 
surface and water vapour in the column. Thus 

P 37 = ( F 37V — T 37H ) 

X exp(0. 1 5 \U + 0.00607 F-4.40), (2.3 ) 

where U is surface wind speed and V is water 
vapour content (SI units), also available from 
SSM/I brightness temperatures. Further details on 
the derivation of the complete and self-consistent 
set of algorithms developed to interpret SSM/I 
data are found in Petty (1990). Similar approaches 
to derive liquid water amount from microwave 
data are reported in Curry etal. (1990) and 
Takeda and Liu (1987). 

It is pertinent to note that the liquid water con- 
tent (LWC) derived from AVHRR as well as from 
microwave data comprise both cloud water and 
precipitating water. Thus there are reasons to 
expect that the analysed amounts should tend to 
be somewhat larger than the vertical integral of the 
simulated water content, which stems only from 
the clouds. We use the notation LWC both for 
analysed and simulated amounts, even if the latter 
ones strictly could be called cloud water content 
(CWC). 


3. Comparison of model clouds and clouds 
derived from satellite measurements 

3.1. Choice of situation 

As mentioned above, there may be difficulties in 
the interpretation of AVHRR data at low sun 
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elevations and over snow covered surfaces. We 
have taken those aspects into account in our 
choice of situation for this first comparison. 
Furthermore, in order to make use of microwave 
data, we have looked for a cloud system that has 
its major parts over the sea surface during the 
period of comparison. This implies that there 
are very few conventional surface observations 
available for verification. We also wanted to find a 
situation with a pronounced non-stationary cloud 
system under motion during the simulation period. 

Thus, we are considering a situation with initial 
time at 26 August 1988 12 GMT and the 30 h 
forecast (till 27 August 1988 18 GMT). The synop- 
tic surface map is shown in Fig. 1. During this 
period, we have AVHRR data centred around 27 
August 1988 8.45 GMT (Fig. 2), which we call 
Savl in the sequel, corresponding to +21 h 
simulation time, and 27 August 1988, 14.30 GMT 
(Fig. 3), Sav2, corresponding to +27 h. We have 
microwave data centred around three different 
points in time, one at 26 August 1988, 19.00 GMT, 


(around +6h), another at 27 August 1988, 
4.30 GMT, (around + 16 h) and a third one 
from 18.30 GMT, (around + 30 h); the last men- 
tioned we call Smi3. It may be pointed out that in 
all those cases, the extracted satellite data cover 
only part of the computational area. 

3.2. Cloud cover from model and AVHRR 
measure men t s - - - - - - 

The cloud cover in the forecast model is defined 
at each level, so we have to calculate a 
2-dimensional field of the cover, corresponding to 
what a satellite sees. There are three basic alter- 
natives for such a calculation, namely, assump- 
tions that clouds at different levels have either 
minimumT maximum or random overlap. There is 
of course some arbitrariness in each of the options. 
According to a recent study by Tian and Curry 
(1989), the two last mentioned alternatives are 
more realistic than the first one. We have adopted 
a maximum/random overlap assumption, imply- 
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Fig. 2. Satellite (NOAA 10) picture in near infrared (channel 2) from passage 27 August 1988 around 8.45 GMT. 
(Satellite Receiving Station, Dundee University, UK.) 


ing that maximum overlap is used as long as the 
cover varies monotonically from level to level in a 
column and that random overlap is taken when 
non-monotonic variation occurs in the same 
column (Sundqvist et al., 1989). 

Figures 4a-b show that there is good overall 
agreement between Savl and the simulated cover 
at the corresponding time. One of the first things 
that we notice is that the satellite analysis has a 
somewhat larger overcast area than the model 
field. This disagreement is seen over the Gulf of 
Bothnia and over the area behind the cold front 
east of Scotland-England where the clear area is 
more extended in the model field than in the 
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analysis. A corresponding deviation is also found 
south-west of Ireland. The synoptic map shows 
that a frontal wave is under development in this 
area. To see those differences quantitatively, we 
have subtracted the model cloud cover from the 
observed one in Fig. 4c. The differences mentioned 
above are dearly seen in this figure. These devia- 
tions, where the analysis shows practically over- 
cast while the model field is cloud- free, could 
possibly be due to an uncertainty in the condensa- 
tion scheme as to what the (threshold) conditions 
should be for condensation to appear. But it seems 
more likely that those deviations are due to the 
model dynamics, which may fail in simulating 
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Fig- 3- Satellite (NOAA 9) picture in near infrared (channel 2) from passage 27 August 1988 around 14.30 GMT. 
(Satellite Receiving Station, Dundee University, UK.) 


realistically the motions connected with (a) the 
fairly weak pressure gradients through a great 
depth o ver ce ntral Sweden and the Gulf of 
’Bothnia, (b) the relatively small scale features 
behind the cold front, and (c) the developing wave 
entering through the southern boundary of the 
""modelarea. 

In the south-west sector of the low, where con- 
vection is taking place (Fig. 2), there are somewhat 
more negative than positive differences (Fig. 4c). 
This is even more evident in the water amount dif- 
fe r enc e as we will see below. The synoptic observa- 
tions in Scotland and Ireland report shower 
activity with cloud cover between three and seven 


octas. These features indicate that the interpreta- 
tion of the AVHRR data has more uncertainty 
over areas with scattered convection. 

Near the northern boundary of the analysis 
region, the difference pattern, with its positive 
values to the north and negative values to the 
south, indicates that the model may have a phase 
error in this condensation system. 

The narrow band of negative values, implying 
more clouds in the model than in the analysis, just 
north and west of the low pressure centre is also 
worth noting in Fig. 4c. Fig. 2 exhibits a narrow 
cloud- free band, which evidently cannot be 
resolved by the model. From Fig. 4b it is seen that 
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the analysis on the model resolution also loses the 
cloud- free feature and shows a partly cloudy band. 

There are several areas with positive values in 
Fig. 4c that are interesting. As examples we may 
particularly regard the area along England’s east 
coast, and the band beginning at 60 °N and west of 
Norway and running toward SSE just west of 
Jutland into Western Germany. In these cases, "the 
model simulates fractional cloudiness ranging 
between 25 and 75%, while the analysis shows a 
cover that systematically is 15-50 units larger. 
These discrepancies are probably due to the 
models parameterization scheme for condensa- 
tion. Thus, it is likely that the scheme tends to 
underestimate the fractional stratiform cloud 
cover. 

We may also make a more statistical com- 
parison by considering the following numbers: in 
the analysis area, the analysed clouds are present 
at 1560 gridpoints with an average cover oT66%, 
while the model clouds are present at 1310 grid- 
points yielding a mean cover of 53 %. Fig. 5 shows 
a frequency distribution of cloud cover categories 
for both the analysed clouds and the model clouds. 


■ Mod-Cover □ Sat-Cover 



CLOUD COVER IN % 

Fig. 5. Frequency distribution of cloud cover categories 
for the analysed clouds and the model clouds. The case is 
That of Sav2 and the forecast is for +27 h. 


We have here chosen the time of Sav2, at which the 
comparison exhibits essentially the same results as 
those shown in Figs. 4a-c. From Fig. 5 we see that 
the analysed and the simulated cloud covers have 
an over all similarity in their bimodal distribu- 
tions. But there are discrepancies, which are most 
pronounced in the extreme categories, that is, for 
zero and for TOO % . There are markedly more 
clbudfree points and fewer overcast points in the 
simulated case than in the satellite analysis. This 
circumstance again shows the model’s tendency to 
underestimate the cloud cover. 

3.3 . Cloud water content from model and A VHRR 
measurements 

The satellite derived liquid water content is com- 
pared with the vertical integral of the predicted 
cloud water content. We noted earlier that the 
2-dimensional cloud cover suffers from some 
arbitrariness due to the necessary assumption 
about overlapping; this is not the case with the 
vertically integrated water content, which is a well 
defined quantity. 

As in the case of cloud cover, there is good 
qualitative agreement between the simulated cloud 
water field and the one obtained from the satellite 
measurements. Figs. 6a-b show that the simulated 
water content has a somewhat wider spread than 
found in Savl, where the higher values appear in 
narrower bands. The maximum values are much 
the same, although there seems to be a few more 
points with the high values (more than 
1.5 kg m” 2 ) in the Savl field (over Western 
Germany). 

Fig. 6c shows the difference between the satellite 
derived water content and the model’s vertically 
integrated cloud water amount. We may first 
observe theTlifference in the convection mentioned 
above in connection with convective cloud cover. 
In Fig. 6c, the simulated water amount is higher 
than the analysed amount over Scotland and 
western Ireland, and over the sea west of there. 
Regarding this discrepancy together with the syn- 
optic reports and cloud cover differences shown 
above, it appears likely that the interpretation of 
the AVHRR data underestimates the convective 
cloud cover and water amount. There is also a 
band along the front from Cornwall toward ENE 
through England where the simulated water 
amount is higher than the analysed one, but in this 
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Fig. 7. The difTerence field, interpreted water amount at the time of Fig. 3, minus simulated cloud water content from 
+ 27hor the model prediction. The scale to the left shows the shading for the different amount intervals unit 
kg-m , (or mm water column). 


■ Cloud Cover □ Cloud Water 
1000 | 



S-Y; M-N Y;Y or N;N S-N; M-Y 
S = Satellite; M = Model; Y = Yes; N = No 


Fig. 8. The figure exhibits the number of grid points of 
agreement and of disagreement between observation and 
simulation for cloud cover and for cloud water content. 


case it is more difficult to judge which of the two 
fields is the more realistic one. 

At the time of Sav2, Fig. 7, there is a more 
distinct difference in the cloud water content, 
with the higher values in the measured field. 
However, there is better agreement in the spatial 
distribution than at the time of Savl. The degree of 
spatial agreement at the time of Sav2 is indicated 
in another way in Fig. 8, which depicts the number 
of grid points of agreement and of disagreement 
between observation and simulation for cloud 
cover as well as for cloud water content. The figure 
also shows that the model underestimates cloud 
cover. Furthermore, the figure indicates that there 
is a discrepancy between the number of points with 
cloud cover and the number of points with cloud 
water. A slight contribution to this discrepancy 
may come from the model output, because the 
fields of cover and water differ by one time step. 
But the major contribution appears to come from 
the analyses, in which we found a distinct dif- 
ference between Savl and Sav2 in spatial spread of 
cloud water. This question is directly connected to 
the matter discussed in the next paragraph. 

There is one feature that is apparent when look- 
ing at Figs. 2 and 3, namely that the clouds appear 
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considerably brighter on the same side of the pic- 
ture as where the sun is, as seen from the satellite. 
This has to do with the anisotropy in the bidirec- 
tional reflectivity, as reported by Taylor and Stowe 
(1984). When calculating the LWC from the 
satellite data, we have regarded the clouds as 
isotropic reflectors, so the anisotropy may explain 
some features that seem somewhat unrealistic. For 
example, from Fig. 6c we see that, in the eastern 
part, the satellite derived LWC at 8.45 GMT is 
higher than the simulated one at 4-21 h. Corre- 
spondingly we see from Fig. 7 that the analysed 
LWC at 14.30 GMT is considerably higher than 
the simulated one at 4-27 h in the western part. 
Those substantial positive deviations may in part 
be explained by the anisotropy, but it is apparent 
that some of them are due to underestimates in 
the model. For instance, the synoptic map at 
09.00 GMT shows a rainband over West Germany 
and the synoptic map at 15.00 GMT shows a rain 
area over southern England and to the south of 
Iceland. The relatively high values of LWC from 
the satellite analysis are thus indirectly 
corroborated by other observations, while the 
simulated precipitation is nearly negligible in the 
same areas. The area of positive difference in 
Fig. 7, extending from about OE, 62N southward 
along the west coast of Norway, may be explained 
by the fact that the precipitating water is not 
included in the model LWC. The simulated 3 h 
precipitation at + 27 h is significant in this area. 
Most of the alternating positive and negative 
differences in the North Sea in Fig. 6c are due to 
the later time of the analysed frontal cloud band 
compared to the simulated cloud band in this area. 

3.4. Cloud water content from model and SSM/f 
measurements 

The LWC derived from the microwave data is 
compared with the vertical integral of predicted 
cloud water content. As indicated in Subsec- 
tion 3.1, we have performed comparisons at three 
different points in time. As the results of those 
comparisons are essentially the same, we have 
chosen to show the results from the 4- 30 h simula- 
tion (Fig. 9a) and Smi3 (Fig. 9b). The difference, 
Smi3 minus the simulation, is exhibited in Fig. 9c. 
The comparison is relevant only over the sea, since 
the microwave interpretation only exists there. It is 
seen that even in the -f-30h forecast, there is 
remarkably good agreement between the positions 


of the frontal bands in Figs. 9a and 9b. An area of 
positive differences over Skagerack west of Sweden 
is seen in Fig. 9c. According to the synoptic map of 
18.00 GMT, it is raining over that area, thus 
indicating t hat th e numerical forecast under- 
estimates the LWC. East of Iceland there is a band 
of positive differences, and a somewhat less 
coherent negative one south of that band. This 
situation indicates that the simulated system has a 
phase lag. However, this may, to some extent, be 
fictitious since the western part of the satellite 
analysis is from the satellite pass at about 
19.30 GMT. 


4. Summary 

In the present study, cloud cover and cloud 
water content from a NWP model have been com- 
pared to the corresponding quantities derived 
from radiances measured by AVHRR and SSM/I 
on satellites. It is demonstrated that this type of 
retrieved meteorological parameters constitute 
potentially useful data samples for quantitative 
verification of model simulations. The results 
obtained in this pilot exercise are indeed encourag- 
ing. But there are some weaknesses in the retrieval 
methods that have to be removed before we can 
have sufficient confidence in the data to make 
quantitative verification meaningful. For instance 
in the retrieval methods applied to AVHRR data, 
it is necessary to treat the anisotropic reflectance 
more accurately. The constraint that the sun eleva- 
tion is above a certain limit should be eliminated. 
The results of the present study indicate that the 
interpreted cloud cover in convective areas may 
have an ^satisfactory uncertainty. Furthermore, 
the inability to distinguish between snow cover 
and clouds is a very important difficulty that has to 
be removed so that the method can be applied with 
confidence at all times of the year. 

As for the liquid water content, there is no 
straightforward answer to what is the best 
estimate, the simulation or the analysis. However, 
by utilizing precipitation reports on synoptic 
maps, we have some possibility to make a judge- 
ment over areas where the simulation and analysis 
have clear differences. 

It is necessary that data from direct 
measurements (e.g., from FIRE and alike) are 
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utilized for calibration of the retrieved cloud 
parameters. To allow calibration of interpretations 
of microwave data, it is desirable that some 
dedicated field experiment is undertaken over 
ocean areas. Cloudiness produced in model 
simulations may also be compared with data from 
in situ measurements; an estimate of the simulated 
water content in the cloud proper is obtained by 
dividing the integrated water amount of a grid 
column with the associated cloud cover. 

Despite uncertainties in the observational data 
considered here, it has yet been possible to infer 
some weaknesses in the model's parameterization 
scheme, by regarding the quantitative differences 
between analysis and simulation. Hence, it is found 
that the scheme has difficulties in handling the 
extremes of cloud cover in stratiform situations; 
there appears to be an underestimate of the cases 
with 100% cover and an overestimate of cloud- 
free areas. 

Extensive utilization of the potential that 
retrieved cloudiness from satellite measurements 
represent will lead to improvements and 
refinements in the interpretation methods. So the 
sooner we begin to make use of those data, the 
sooner we will acquire a satisfactory level of 
accuracy and confidence in the retrieved data. 
In turn, this will facilitate improvements and 
refinements in the condensation-cloud para- 
meterization schemes. As our parameterization 
schemes improve, we can utilize them to derive the 
vertical distribution of condensation and cloud 
processes. It will then be possible to improve the 
humidity analysis as well as to establish initial 
states of release of latent heat and of cloud water 
content. A confident description of the latent 
heating is of great value in the initialization 
process. As more advanced condensation-cloud 
schemes become adopted in prediction models, it 
will be of great importance (e.g., for minimizing 
spin-up) that a realistic initial state of cloud water 
content is given. 
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6. Appendix 

List of ACRONYMS 

AMUB Section of Meteorology, University of 

Bergen 

AVHRR Advanced Very High Resolution 

Radiometer 

BSC Bergen Scientific Centre (IBM ) 

CFL Courant-Friedrichs-Levy (numerical 

stability condition) 

CWC Cloud Water Content 

FIRE First ISCCP Regional Experiment 

GCM General Circulation Model 

• ISCCP International Satellite Cloud Climatol- 

ogy Project 

LWC Liquid Water Content 

NRSC Nansen Remote Sensing Center 

NWP Numerical Weather Prediction 

SMMR Scanning Multichannel Microwave 

Radiometer 

SSM/I Special Sensor Microwave Imager 
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